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PREFACE 


This  final  report  describes  the  work  performed  by  Technion  Research 
and  Development  Foundation  Ltd.,  Technion,  Haifa  32000,  Israel,  under 
U.S.  Air  Force  Contract  F49620-83-C-0057.  The  report  covers  the  period 
from  January  1983  through  August  1985. 

This  study  of  noncontacting  coned-face  mechanical  seal  dynamics  was 
sponsored  by  the  Aero  Propulsion  Laboratory  of  the  Air  Force  Wright 
Aeronautical  Laboratories  (AFWAL),  Wright-Patterson  Air  Force  Base,  Ohio 
45433-6563  under  project  3048,  "Fuels,  Lubrication  and  Fire  Protection", 
Task  304806,  "Aerospace  Lubrication",  Work  Unit  30480629,  “Mechanical 
Face  Seal  Dynamics".  Lt  G.  Gainer  and  R.  Dayton  were  the  USAF 
AFWAL/POSL  Project  Engineers.  The  work  at  Technion  was  performed  by 
Dr  I.  Etsion,  who  was  the  Principal  Investigator,  and  Dr  1.  Green. 
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SECTION  I 


INTRODUCTION 

Radial  face  seals  are  used  in  a  host  of  applications  to  seal  liquids  of 
various  types.  These  liquids  range  from  lubricants  to  highly  toxic  chemicals 
and  acids.  The  applications  range  from  helicopter  transmissions  to  nuclear 
reactor  cooling  pumps  and  submarine  propeller  shafts.  The  function  of  these 
seals  is  to  restrict  the  leakage  of  the  sealed  fluid  and  to  prevent  the  entry 
of  solid  and  liquid  debris.  A  seal  permits  a  rotating  shaft  to  penetrate  an 
enclosure  (transmission  box,  submarine  hull,  pump  housing,  etc.)  while  main¬ 
taining  separation  of  the  environments  on  the  inside  and  outside  of  the  enclo 
sure.  To  do  this,  one  component  of  the  seal  is  attached  to  the  shaft  and 
rotates  with  it,  and  the  other  component  is  attached  to  the  housing  and  is 
non-rotating.  One  of  the  two  components  is  flexibly  mounted  to  provide  angu¬ 
lar  and  axial  freedom  of  motion.  Mechanical  forces  and  fluid  film  pressures 
should  tend  to  force  the  flexibly  mounted  face  into  alignment  with  the  other 
component.  Relative  sliding  motion  takes  place  between  the  faces  of  these  two 
components.  In  order  to  avoid  wear  and  to  achieve  long  life,  these  two 
sealing  faces  must  be  separated  by  a  film  of  the  sealed  liquid.  This  film 
must  be  very  thin  to  keep  the  leakage  rate  within  acceptable  limits.  Thus, 
the  requirements  of  lubrication  and  leakage  tend  to  be  conflicting. 

There  have  been  many  hypotheses  put  forth  to  explain  the  mechanisms 
responsible  for  the  development  of  the  lubricating  film  pressure  that  acts  to 
separate  the  primary  seal  faces.  These  hypotheses  include  the  following: 
surface  angular  misalignment,  surface  waviness,  surface  asperities,  vaporiza¬ 
tion  of  the  fluid  film,  axial  vibration,  and  thermal  deformation.  A  good 


review  of  the  state-of-the-art  is  given  in  Reference  1,  and  some  of  the  latest 
results  on  seal  theory  are  presented  in  References  2-8.  Examination  reveals  that 
these  published  theories  are  of  very  limited  use  from  an  operation  prediction 
standpoint.  Briefly  stated,  face  seal  lubrication  theory  is  very  primitive  as 
compared  to  journal  bearing  theory.  Classical  journal  bearing  theory  applies 
to  face  seal  lubrication,  but  seal  dynamics,  which  is  thought  to  be  of  major 
importance,  is  poorly  understood. 

Dynamic  instability  in  the  form  of  vibration  of  flexibly  mounted  elements 
was  experimentally  observed  in  face  type  seals  (References  9-13).  However, 
only  a  few  attempts  have  been  made  so  far  to  analyze  seal  dynamics.  Most  of 
these  analyses  overlook  some  of  the  fluid  film  or  flexible  support  effects 
(References  14-18).  A  fuller  treatment  of  the  problem  is  presented  in 
Reference  19,  but  it  treats  a  special  case  and  does  not  offer  general  design 
criteria.  In  more  recent  work  (References  20-22)  the  motion  of  a  flexibly 
mounted  ring  of  a  noncontacting  face  seal  is  described  in  its  three  major 
degrees  of  freedom  (one  axial  and  two  angular).  Hydrodynamic,  hydrostatic,  and 
squeeze  film  effects,  as  well  as  the  contribution  of  the  springs  of  the  flexible 
support,  are  considered  and  it  is  shown  how  the  seal  stability  is  affected  by 
various  design  parameters.  The  analysis  in  References  20-22  is,  however, 
somewhat  limited  in  that  it  does  not  contain  the  effect  of  an  elastomeric  secon¬ 
dary  seal  which  often  exists  in  face  seals.  This  element  adds  both  stiffness 
and  damping  to  the  flexible  support  but  the  stiffness  and  damping  coefficients 
are  not  constants.  In  the  case  of  seal-rotor  runout,  the  stiffness  and  damping 
coefficients  of  the  elastomer  are  frequency  dependent  which  complicates  the 
solution  of  the  dynamic  problem. 

The  general  thrust  of  this  work  is  aimed  at  studying  the  dynamic  properties 
of  elastomeric  secondary  seals  and  incorporating  their  effects  in  a  dynamic  ana¬ 
lysis  of  mechanical  face  seals. 


SECTION  II 


APPROACH 


The  effort  was  divided  into  three  main  tasks.  Task  1  was  devoted  to 
studying  experimentally  the  dynamic  properties  of  elastomeric  0-rings.  A  test 
rig  to  simulate  the  behavior  of  secondary  seal  0-rings  in  a  practical  mechanical 
seal  was  built.  Stiffness  and  damping  coefficients  of  0-rings  of  various  sizes 
and  materials  were  measured  for  a  range  of  frequencies  and  amplitudes.  This 
effort  is  described  in  Appendix  A. 

Task  2  was  devoted  to  developing  the  kinematic  model  of  a  mechanical  seal 
and  studying  the  effect  of  anti-rotation  locks  (Appendix  B).  Following  this 
study  the  dynamic  behavior  of  a  non-contacting,  coned-face  mechanical  seal  was 
analyzed.  The  analysis  was  based  on  linearization  of  the  equations  of  motion  of 
the  flexibly  mounted  element  of  the  seal.  This  approach  yielded  analytical 
expressions  and  enabled  parametric  investigation  of  the  stability  threshold  and 
steady  state  response  of  the  seal  (see  Appendix  C). 

During  the  third  task,  the  complete  non-linear  dynamic  analysis  of  the  seal 
was  performed.  A  numerical  code  was  developed  to  solve  the  non-linear  equations 
of  motion  of  the  flexibly  mounted  element.  A  comparison  was  made  between  the 
results  obtained  from  the  non-linear  analysis  and  those  obtained  analytically 
from  the  small  perturbation  analysis  of  Task  2.  This  effort  is  described  in 
Appendix  D. 


SECTION  III 


MAIN  RESULTS  AND  CONCLUSIONS 


1.  In  practical  seals,  the  dynamic  properties  of  the  secondary  seals'  0-rings 
are  negligibly  affected  by  seal  runout.  Increasing  the  shaft  speed 
results  in  higher  stiffness  coefficient  of  elastomer  0-rings  but  may  either 
increase  or  decrease  the  damping  coefficient  depending  on  the  elastomer 
material.  The  stiffness  K  and  damping  D  can  be  expressed  by  an  exponen¬ 
tial  function 

K  =  Atj® 

D  =  aoj^ 

where  w  is  the  shaft  speed  and  A,  B,  a,  and  b  are  constants  related  to 
elastomer  material  and  geometry.  Stiffness  provided  by  the  secondary  seal 
0-rings  is  generally  much  higher  than  that  of  mechanical  springs. 

2.  A  stability  threshold  can  be  found  in  terms  of  a  critical  shaft  speed  below 
which  stable  seal  operation  is  assured.  This  critical  speed  is  normally 
very  high  and,  hence,  stability  threshold  should  not  be  a  problem  in  most 
practical  applications.  It  may  become  a  problem  in  cases  where  the  shaft 
speed  is  very  high,  stator  has  a  large  mass  and  sealed  pressure  is  low. 

3.  Non-contacting  seals  operate  with  relative  misalignment  between  the  mating 
faces.  This  relative  misalignment  depends,  among  other  parameters,  on  the 
runout  of  the  rotor  and  is  usually  a  function  of  time  too,  varying  cycly- 
cally  at  the  shaft  frequency. 
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A  critical  rotor  runout  exists  which,  if  exceeded,  may  cause  failure  by 
local  rubbing  contact  between  the  mating  faces  of  the  seal  due  to  a  too  high 
relative  misalignment. 

Failure  because  of  an  excessive  leakage  resulting  from  an  increase  in  seal 
clearance  due  to  rotor  runout  is  likely  to  occur  much  before  the  critical 
runout  is  reached.  This  can  happen  in  seals  with  low  stiffness  of  the 
fluid  film  and  flexible  support.  Small  design  clearance,  high  sealed 
pressure  and  optimum  coning  provide  high  stiffness  of  the  fluid  film  and 
reduce  the  danger  of  a  too  high  operating  clearance. 

A  comparison  of  the  results  obtained  from  a  small  perturbation  analysis 
and  a  full  non-linear  analysis  shows  a  very  good  correlation  for  most 
cases  of  practical  applications.  Fair  correlation  exists  even  in  cases 
which  are  clearly  out  of  the  range  of  small  perturbation.  In  these 
cases,  the  more  simple  small  perturbation  analysis  yields  results  that 
are  on  the  safe  side  with  regard  to  the  critical  rotor  runout.  Hence,  the 
analytical  results  of  the  small  perturbation  analysis  can  be  used  safely 
over  the  full  range  of  design  parameters  and  operation  conditions. 
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STIFFNESS  AND  DAMPING  CHARACTERISTICS  OF  ELASTOMER  0-RINGS 
SECONDARY  SEALS  SUBJECTED  TO  RECIPROCATING  TWIST 

I.  Green  and  I.  Etsion 
Department  of  Mechanical  Engineering 
Technion,  Israel  Institute  of  Technology 
Haifa,  Israel 


Summary 

A  test  is  described  which  simulates  forced  vibration  conditions  of  an  0-ring 
secondary  seal  in  a  mechanical  face  seal.  The  0-ring  cross  section  is  subjected  to  a 
reciprocating  harmonic  twist  similar  to  that  resulting  from  a  rotor  runout  in  an  actual 
seal.  Stiffness  and  damping  coefficients  of  various  elastomer  0-rings  were  measured  and 
are  presented  for  a  range  of  forcing  amplitudes  at  frequencies  from  100  to  300  Hz.  The 
effects  of  elastomer  material,  0-ring  geometry,  and  vibration  amplitude  and  frequency 
are  discussed. 
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=  constants 

=  damping  coefficient  ,  Ns/m 
=  force ,  N 

=  stiffness  coefficient,  N/m 
=  complex  stiffness  of  elastomer 
=  mass  ,  kg 
=  time,  sec 

=  mass  displacement  relative  to  base,  m 
=  mass  amplitude  relative  to  base,  m 
=  base  displacement,  m 
=  base  amplitude,  m 
=  transmissi bi lity , 

=  loss  coefficient,  Ow/K 
=  phase  shift,  rad 
=  frequency,  rad/s 
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1.  INTRODUCTION 

A  mechanical  face  seal  (Fig.  1)  is  a  complex  dynamic  system.  The  rotating 
seat  with  its  unavoidable  axial  runout  forces  both  axial  and  angular  vibrations  of  the 
flexibly  mounted  ring.  This  can  be  schematically  described  by  the  simplified  single 
degree  of  freedom  model  in  Fig.  2.  The  harmonic  motion  of  the  seat  is  transmitted  via 
the  fluid  film  to  the  flexibly  mounted  ring  represented  by  the  mass  m.  The  flexible 
support  of  the  seal  ring  consists  of  springs  in  conjunction  with  a  secondary  seal  ele¬ 
ment.  In  many  cases  this  secondary  seal  is  an  elastomer  0-ring  having  stiffness  K  and 
damping  D  which  affect  the  dynamic  behavior  of  the  entire  seal  system.  In  Ref.  ®(1) 
an  order  of  magnitude  of  the  stiffness  and  damping  coefficients  of  various  components 
of  the  seal  system  is  assessed.  The  0-ring  stiffness  is  shown  to  be  about  two  orders  of 
magnitude  less  than  the  fluid  film  stiffness  but  about  an  order  of  magnitude  higher 
than  the  springs  stiffness.  The  0-ring  damping  is  estimated  about  four  orders  of  magni¬ 
tude  less  than  the  fluid  film  damping,  but  this  is  the  only  damping  provided  by  the 
flexible  support. 

Inspite  its  relative  importance,  the  effect  of  the  elastomer  0-ring  on  face 
seal  dynamics  was  not  considered  in  seal  analyses  published  so  far.  A  step  in  the  right 
direction  was  made  in  Ref.  (2).  Axial  forces  transmitted  by  0-rings  subjected  to  a  re¬ 
ciprocating  drag  were  measured  for  various  amplitudes  and  frequencies.  However,  for  a 
general  dynamic  analysis  the  stiffness  and  damping  coefficients  of  the  0-ring,  not  the 
forces,  are  required. 

An  excellent  source  on  dynamic  properties  of  elastomers  is  provided  in 
Ref.  (3).  A  main  conclusion  emerging  from  this  source  is  that  no  solid  theoretical  pre¬ 
diction  of  elastomers  dynamic  characteristics  is  available  at  present.  Stiffness  and 
damping  coefficients  can  at  best  be  measured  and  are  found  to  be  affected  by  composi¬ 
tion,  geometry,  frequency,  strain,  temperature,  and  preload. 

While  a  great  deal  of  information  has  been  gathered  on  various  elastomer 
springs  and  dampers  this  is  not  the  case  for  0-rings.  Refs.  (4)  and  (5)  present 
results  on  measured  dynamic  properties  of  elastomer  0-rings  under  reciprocating  radial 
squeeze,  similar  to  the  mode  of  application  in  squeeze  film  dampers.  The  secondary  seal 
0-ring  in  a  mechanical  face  seal  is,  however,  subjected  to  a  reciprocating  twist  of  its 
cross  section.  This  mode  of  operation  is  the  result  of  the  flexibly  mounted  ring  track¬ 
ing  the  axial  runout  of  the  rotating  seat. No  information  on  the  dynamic  characteristics 
of  0-rings  under  reciprocating  harmonic  twist  of  their  cross  section  is  available  in 
the  literature.  This  report  describes  an  effort  to  obtain  such  missing  information 
which  is  vital  for  better  understanding  of  mechanical  seals  dynamics. 


2.  BACKGROUND 

Elastomers  are  characterized  by  stress  relaxation  whereby  the  load  required 
to  maintain  a  constant  strain  decreases  with  time.  If  a  load  cycle  is  applied  repeated¬ 
ly  to  an  elastomer  element,  the  process  of  relaxation  causes  hysteresis  which  in  turn 
results  in  loss  of  energy  per  each  cycle.  Hence,  elastomers  under  cyclic  loading  pro¬ 
vide  energy  dissipation  or  hysteretic  damping. 

The  force  and  displacement  in  a  linear  elastomer  element  undergoing  a  repea¬ 
ted  deflection  cycle  x  =  x*ei‘^t  may  be  related  by  (Ref.  3). 

F  =  k*x  (1) 

where  k*  is  a  complex  number  in  the  form 

k*  =  ki  +  ikj,  (2) 

and  is  generally  a  function  of  frequency,  temperature,  dynamic  strain  and  geometry,  as 
well  as  material  composition.  The  quantities  ki  and  kj  are  generally  referred  to  as  the 
stiffness  and  damping  of  the  elastomer.  The  damping  ability  of  an  elastomer  element  is 
sometimes  expressed  in  terms  of  a  loss  coefficient  n  which  is  a  nondimensional  ratio  of 
the  energy  dissipated  per  cycle  to  the  energy  stored  at  maximum  displacement 


Various  test  methods  for  determining  the  dynamic  properties  of  elastomers  are 
described  in  Ref.  p).  Of  all  these  methods  the  most  advanced  one  is  the  Base  Excita¬ 
tion  Resonent  Mass  (BERM)  test  method,  which  is  described  schematically  in  Fig.  3.  By 
measuring  independently  the  motion  of  the  base  y  =  y{t),  the  motion  of  mass  x  =  xU) 
relative  to  the  base,  and  the  phase  shift  (>  between  these  two  motions,  the  stiffness 
and  damping  of  the  elastomer  can  be  found.  Assuming  the  Kelvin-Voigt  model  for  the 
elastomer  as  shown  in  Fig.  3,  we  have 

m(x+y)  +  Dx  +  Kx  =  0  (4) 

where  D  and  K  are  the  damping  and  stiffness  coefficients,  respectively,  of  the  elasto¬ 
mer.  Assuming  a  base  excitation  of  the  form 

y  =  y^sinCwt+o)  (5) 

and  a  mass  response  relative  to  the  base 


X  =  x^sinut 


eq.  (4)  becomes 

-mw^x^si  nost  +  DwXjjCOSwt  +  Kx^sinwt  =  mw^y^si  n((i)t+(}))  (7) 

Equating  terms  containing  sinwt  and  coswt  separately  we  find 

D  =  tiv^  (8) 

K  =  noj^  +1)  (9) 

where  a  is  the  ratio  of  the  relative  response  amplitude  x  to  the  base  excitation  am¬ 
plitude  Vq,  is  the  phase  shift  between  x  and  y  and  u  is°the  frequency  of  excitation. 

The  stiffness  and  damping  coefficients  K  and  D  are  related  to  ki  and  kz  of 
eq.  (2)  in  the  form 

ki  =  K  (10a) 

kj  =  Du  (10b) 

Hence,  by  eqs.  (8)  to  (10),  the  loss  coefficient  n  of  eq.  (3)  is 

n  =  sincj) 

'  bi+ CO  S<{) 

3.  EXPERIMENTAL  SETUP  AND  TEST  PROCEDURE 

A  test  rig  was  set  up  as  shown  in  Fig.  4  to  obtain  the  stiffness  and  damping 
coefficients  of  various  elastomer  0-rings  subjected  to  reciprocating  twist  of  their 
cross  section.  A  holder  (1)  attached  to  the  vibrating  table  (2)  of  an  electromagnetic 
shaker  is  holding  the  base  (3)  by  means  of  the  screw  (4).  Three  different  base  sizes 
were  used  to  accomodate  an  0-ring  (5)  of  76.2  mm  mean  outer  diameter  with  three  diffe¬ 
rent  nominal  cross  section  diameters  of  1.78  im,  3.53  mm,  and  5.33  mm  (1/16".  1/8", 
and  3/16").  A  mass  of  1  kg  in  the  form  of  an  annular  ring  (6)  was  fitted  on  the  0-ring. 
The  size  of  the  0-ring  groove,  and  the  inner  diameter  of  the  ring  were  selected  from 
ref.  (6)  for  industrial  0-ring  static  seals.  This  provided  a  31  percent  squeeze  for 
the  1/16"  0-ring  and  20  percent  squeeze  for  both  the  1/8"  and  3/16"  0-rings.  Two  M60 
Dymac  eddy  current  proximity  probes  (7)  and  (7a)  were  used  to  measure  the  base  motion 
y=y(t)  and  the  mass  motion  relative  to  the  base  x  =  x(t),  respectively. 

A  date  acquisition  system  based  on  a  PDP  11/40  minicomputer  was  used  to  re¬ 
cord  the  probes  output  and  to  analyze  the  data.  The  two  signals  y(t)  and  x(t)  were 


sampled  at  a  rate  of  1000  Hz  after  proper  filtration  and  amplification.  Tests  were 
conducted  at  three  different  amplitudes,  y  ,  of  the  vibrating  base.  The  amplitudes 
were  13  ym,  16  ym,  and  20  ym.  These  amplitudes  are  small  enough  to  prevent  slippage 
of  the  mass  ring  on  the  0-ring  but  are  of  the  same  order  of  seal  seat  motion  due  to 
seat  runout  in  practical  face  seals.  It  was  important  to  prevent  slippage  of  the  mass 
ring  in  order  to  eliminate  the  effect  of  coulomb  friction  on  the  dynamic  properties  of 
the  0-ri  ngs . 

Ten  different  frequencies  in  the  range  100  Hz  to  300  Hz  were  selected  for  the 
tests.  At  each  one  of  the  preselected  frequencies  the  shaker  amplitude  y  was  main¬ 
tained  constant.  The  output  of  the  probe  monitoring  the  base  motion  was  “observed  on  a 
digital  voltmeter  to  facilitate  control  of  the  shaker  amplitude.  Output  of  the  two 
proximity  probes  was  also  observed  on  an  asciloscop  to  detect  any  possible  unusual  be¬ 
havior  of  the  test  rig  (like  mass  slippage,  for  example).  Only  test  points  which 
looked  satisfactory  were  recorded  and  stored  in  thecomputer.  This  procedure  was  repea¬ 
ted  for  each  of  the  three  different  amplitudes  y  .  Several  tests  were  rerun  after 
taking  the  0-ring  off  the  rig  and  reassembling  again  in  order  to  check  the  repeatabili¬ 
ty  of  the  results.  All  tests  were  run  at  room  temperature  which  was  about  25°  C. 

A  spectral  analysis  code  was  used  to  find  the  discrete  frequency  w  at  each 
test  point  along  with  the  corresponding  transmissi bi lity  a,  and  phase  shift  (|i.  Equa¬ 
tions  (8),  (9)  and  (11)  were  then  used  to  calculate  the  stiffness,  damping,  and  loss 
coefficients,  respectively.  As  a  final  step  a  least  square  procedure  was  utilized  to 
find  a  best  fit  for  the  collection  of  data  points.  This  provided  expressions  for  the 
stiffness  and  damping  as  functions  of  the  frequency,  at  a  constant  amplitude  y-j-in  the 
exponential  form 


K  =  Alj® 

(12) 

n 

0  -  aoj 

(13) 

where  A,  B,  a,  and  b  are  constants  related  to  the  elastomer  material  and  to  the  0-ring 
geometry.  The  expression  for  the  loss  coefficient,  according  to  eq.  (11)  and  eqs.  (12), 
(13),  is  simply 

-  a  ,  (b-B+1) 

4.  RESULTS  AND  DISCUSSION 

Results  of  the  stiffness,  damping  and  loss  coefficients  of  Nitrile  (Buna  N) 
and  Fluorocarbon  (Viton  75)  elastomer  O-rings  are  presented  in  Figs.  5,  6  and  7  re¬ 
spectively.  Varying  the  amplitude  yg  in  the  investigated  range  which  was  IS^msy^  iZOum 
changed  the  dynamic  properties  by  less  than  Z%.  Hence,  no  distinction  is  made  in  Figs. 

5  to  7  between  the  various  amplitudes.  Scatter  of  the  test  results  after  reassembling 
the  test  rig,  as  described  in  the  previous  section,  was  up  to  30X  as  can  be  seen  from 
Figs.  5a  and  6a  in  the  case  of  the  Buna  N  with  1/8"  cross  section  diameter.  It  should 
be  noticed,  however,  that  for  each  given  0-ring  assembly  the  behaviour  was  consistent 
and  the  data  points  were  lined  quite  accurately  along  a  straight  line  on  the  log  -  log 
scale,  obeying  the  relations  (12)  and  (13).  Apparently  the  stiffness  and  damping  are 
sensitive  to  some  factor  related  to  the  mode  of  assembling  0-ring  in  the  test  rig 
(which  could  be  initial  twist,  for  example). 

The  effect  of  0-ring  cross  section  diameter  is  the  same  in  both  materials. 

The  smallest  diameter  provides  the  highest  stiffness  and  damping.  This  is  partly  attri¬ 

buted  to  the  higher  squeeze  which  was  31  percent  for  the  1/16"  diameter  0-ring  compared 
to  20  percent  for  the  1/8"  and  3/16"  diameters.  However,  most  of  the  difference  in 
stiffness  and  damping  is  due  to  the  change  in  diameter  as  can  be  seen  from  the  results 
of  the  two  larger  diameter  0-rings. 

The  effect  of  elastomer  material  is  very  distinct.  In  general  the  Viton  75 
has  stiffness  that  is  about  twice  the  stiffness  of  Buna  N  for  the  same  frequency 

(Fig.  5).  The  damping  of  Viton  75  is  about  five  times  larger  than  that  of  Buna  N  in  the 

frequency  range  of  the  tests  (Fig.  6). 

Increasing  the  frequency  increases  the  stiffness  in  both  Viton  75  and  Buna  N 
(Fig,  5).  The  damping  coefficient,  however,  behaves  differently  for  the  two  materials 


and  while  an  increase  in  frequency  increases  the  damping  of  Buna  N  (Fig.  6a)  it  redu¬ 
ces  the  damping  of  Viton  75  (Fig.  6b).  A  reduction  in  damping  coefficient  with  increa¬ 
sing  frequency  was  also  observed  in  Refs.  (4)  and  (5)  for  various  elastomer  materials 
when  0-rings  were  subjected  to  radial  squeeze  oscilation.  The  loss  coefficient  of  the 
two  materials  increases  with  the  frequency  (Fig.  7). 

Table  1  sunmarizes  the  results  in  terms  of  the  constants  A,  B,  a,  and  b  of 
equations  (12)  and  (13).  In  general,  for  the  range  of  frequencies  tested,  the  stiffness 
coefficient  of  Viton  75  was  between  0.7  MfJ/m  to  4  fWm  and  its  damping  between  0.8  to 
5  kN  s/m.  The  stiffness  coefficient  of  the  Buna  N  was  between  0.4  to  2  MN/m  and  its 
damping  coefficient  between  0.25  to  0.8  kN  s/m.  These  values  are  of  the  same  order  of 
magnitude  that  was  assessed  in  Ref.  (1).  The  stiffness  and  damping  coefficients  repor¬ 
ted  here  are  also  of  the  same  order  of  magnitude  of  these  reported  in  Ref.  (4)  for 
Viton  70  and  Buna  N  0-rings  of  similar  geometry  under  reciprocating  radial  squeeze. 

5.  CONCLUDING  REMARKS 

A  test  program  for  measuring  stiffness  and  damping  coefficients  of  elastomer 
0-rings  is  described.  The  0-rings  are  subjected  to  a  reciprocating  twist  of  their  cross 
section,  similar  to  the  mode  of  operation  of  secondary  seals  in  mechanical  face  seals. 
Effects  of  elastomer  material,  geometry,  frequency  and  amplitude  of  vibration  were 
examined.  It  was  found  that  in  the  range  of  tested  frequencies  and  base  amplitudes 
(100  to  300  Hz,  and  13  to  20  pm)  the  amplitude  effect  on  the  dynamic  characteristics 
is  negligible.  In  practice  the  amplitude  depends  on  the  amount  of  runout.  Hence,  it 
may  be  concluded  that  in  practical  seals  the  dynamic  properties  of  the  secondary 
seals  0-rings  are  independent  of  the  amount  of  seal  runout. 

Increasing  the  frequency  of  vibration  results  in  a  higher  stiffness  coeffici¬ 
ent  of  elastomer  0-rings  but  may  either  increase  or  decrease  the  damping  coefficient 
depending  on  the  elastomer  material.  Larger  cross  section  diameters  provide  lowerstiff¬ 
ness  and  damping. 

Both  stiffness  and  damping  coefficients  can  be  expressed  by  an  exponential 

function 


K  =  Atu® 


where  A,  B,  a  and  b  are  constants  related  to  elastomer  material  and  geometry.  In  general 
the  stiffness  provided  by  a  secondary  seal  0-ring  may  be  much  higher  than  that  of  the 
mechanical  springs  found  in  the  flexible  support  of  mechanical  seals.  The  damping  pro¬ 
vided  by  the  0-ring  is  the  only  damping  in  the  support.  Hence,  dynamic  properties  of 
secondary  seals  0-rings  have  to  be  thoroughly  considered  in  seal  analysis  and  design. 
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Table  1  summary  of  the  results  for  Buna  N  and  Viton  75  0-rings. 
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ABSTRACT 


A  kinematic  model  of  mechanical  face  seals  is  presented.  Two  basic  seal 
arrangements  are  considered;  A  flexibly  mounted  stator  with  antirotation 
locks,  and  a  flexibly  mounted  rotor  with  positive  drive  devices.  The  equation 
of  kinematic  constraint  is  derived  and  presented  in  a  simple  form  for  all  the 
possible  types  of  antirotation  or  positive  drive  mechanisms  found  in  practical 
seals.  This  simple  form  is  then  used  to  derive  the  dynamic  moments  acting  on 
the  flexibly  mounted  element  of  the  seal. 


NOMENCLATURE 


I  transverse  motnent  of  inertia 

1 2  polar  moment  of  inertia 

L  relative  angular  momentum 

T  transmission  law 

T  dynamic  moment 

a  =  -111 
Y  nutation 

X  absolute  angular  velocity 

ill  relative  precession 

l|lr  absolute  rotor  precession 

ijig  absolute  stator  precession 

relative  spin 

ta)  shaft  angular  velocity 

(l  relative  angular  velocity 

cDc  reference  angular  velocity 


Subscripts 


INTRODUCTION 


Seal  dynamics  has  become  the  subject  of  many  investigations  in  the  last 
decade  [l] .  Much  effort  is  devoted  to  analyze  the  time-dependent  behavior  of 
the  flexibly  mounted  element  of  the  seal.  This  element  can  be  either  the 
rotating  one,  as  in  many  low  speed  applications,  or  the  stationary  one  as  shown 
in  Fig.  1.  Its  motion  is  affected  by  factors  such  as  axial  runout,  shaft 
vibration,  dynamic  properties  of  both  the  flexible  support  and  the  lubricating 
fluid  film,  etc.  These  factors  were  considered  in  previous  works  [l,3,4]  and 
the  existing  theoretical  models  are  quite  close  to  the  realistic  seals.  There 
is,  however,  one  aspect  which  was  overlooked  so  far,  but  nevertheless  plays  an 
important  role  in  seal  dynamics.  This  is  the  constraint  imposed  on  the 
flexibly  mounted  element  by  the  positive  drive  devices  in  the  case  of  a 
flexibly  mounted  rotor  or  by  the  anti-rotation  locks  (see  Fig.  1)  in  the  case 
of  a  flexibly  mounted  stator.  Understanding  of  this  constraint  is  essential 
for  a  correct  formulation  of  the  seal  kinematics.  Unfortunately,  there  are 
numerous  different  arrangements  of  positive  drive  or  anti-rotation  devices 
[5,6],  e.g.  dents,  keys,  pins,  slots  and  ears,  and  bellows  to  name  just  a  few. 
In  addition,  the  number  of  units  in  a  particular  arrangement  may  vary  in  dif¬ 
ferent  designs,  ranging,  for  example,  from  one  to  four  pins  per  seal. 
Manufacturing  tolerances  regarding  these  devices  are  fairly  large  and  there¬ 
fore  even  in  cases  where  several  drives  or  locks  are  present,  only  one  of  them 
may  actually  be  effective. 


The  constraint  situation  described  above  complicates  any  attempt  to  deal 
accurately  with  the  kinematic  model  of  mechanical  seals.  Such  a  model  is, 
however,  necessary  for  the  derivation  of  the  equations  of  motion.  The  present 
paper  describes  a  general  treatment  that  offers  a  fairly  accurate  solution  to 
this  complex  problem.  Based  on  the  fact  that  the  angular  displacements  of  the 
flexibly  mounted  element  are  very  small,  it  will  be  shown  that  a  first  order 
approximation  serves  as  a  good  general  model,  with  a  truncation  error  of  order 
where  y  «  1.  Finally,  the  dynamic  moments  that  act  upon  the  flexibly 
mounted  element  will  be  derived  for  the  two  basic  arrangements  where  this  ele¬ 
ment  is  either  the  rotor  or  the  stator. 


THE  KINEMATIC  MODEL 


Fig.  2  presents  schematically  a  model  of  a  mechanical  face  seal  which 
assists  in  understanding  the  kinematics  of  the  flexibly  mounted  seal  element. 

The  figure  shows  a  system  of  two  rings.  An  outer  ring  to  which  a  reference  XYZ 
is  attached,  and  an  inner  ring  with  two  slots  to  accomodate  two  pins  that  are 
fixed  to  the  outer  ring  along  the  Y  axis.  The  inner  ring  represents  the 
flexibly  mounted  element  of  the  seal  and  is  free  to  have  two  orthogonal  tilts 
about  two  of  its  diameters. 

Two  cases  will  be  considered.  1)  The  flexibly  mounted  seal  element  is  sta¬ 
tionary,  and  2)  the  flexibly  mounted  seal  element  is  rotating.  In  the  first 
case  the  outer  ring  represents  the  seal  housing  (see  Fig.  1),  the  two  pins 
represent  the  antirotation  locks,  and  the  reference  XYZ  is  inertial.  In  the 
second  case  the  outer  ring  represents  the  shaft,  the  two  pins  represent  the 
positive  drive  mechanism  and  the  reference  XYZ,  together  with  the  outer  ring, 
rotates  at  an  angular  velocity  ui  about  the  Z  axis. 

The  resultant  of  the  two  tilts  of  the  inner  ring  can  be  described  by  the 
two  Eulerian  angles  y  and  Ijl  (see  Fig.  2).  The  angle  y  is  the  nutation  of  the 
inner  ring  about  the  axis  x  of  a  reference  system  xyz.  This  reference  system 
is  free  to  rotate  with  respect  to  the  inner  ring  so  that  the  axis  y  is  always 
directed  to  the  point  of  maximum  distance  from  the  XY  plan.  The  angle  l|l  is  the 
precession  of  the  reference  xyz  with  respect  to  the  reference  XYZ.  The  axis  z 
of  the  rotating  xyz  reference  coincides  with  the  principal  axis  of  the  inner 


ring.  It  is  this  axis  about  which  the  inner  ring  has  a  spin  (|)  with  respect  to 
the  xyz  reference. 

An  observer  located  in  the  reference  xyz  sees  the  reference  XYZ  and,  hence, 
the  outer  ring  rotating  through  an  angle  -Ijl  about  axis  Z  while  the  inner  ring 
rotates  through  an  angle  (|)  about  axis  z.  The  kinematic  constraint  forces  the 
two  rings  to  complete  one  revolution  simultaneously,  while  any  pair  of 
corresponding  points  on  the  circumference  of  the  two  rings  return  to  their  ori¬ 
ginal  relative  positions  after  the  completion  of  each  revolution.  This  kinema¬ 
tic  quality  is  characteristic  of  any  universal  joint  and,  hence,  the  seal  model 
of  Fig.  2  can  be  represented  by  a  universal  joint  as  shown  in  Fig.  3.  Here,  the 
rotation  6  =  -ill  is  the  input  to  the  joint  related  to  the  outer  ring,  and  the 
rotation  (|)  is  the  output  from  the  joint  related  to  the  inner  ring. 

The  kinematic  constraint  represented  by  the  two  pins  in  Fig.  2  reduces  the 
number  of  rotational  degrees  of  freedom  of  the  system  into  two,  and  dictates  a 
certain  relation  between  the  Eulerian  angles.  This  relation,  known  as  the 
equation  of  kinematic  constraint,  has  the  general  form 

^  (1) 

and  is  typical  for  universal  joints  (see  for  example  the  case  of  a  Hooke  Joint 
in  [yj  pp.  270-272). 

As  stated  in  the  Introduction,  numerous  different  arrangements  of  antirota¬ 
tion  locks  and  positive  drive  mechanisms  can  be  found  in  mechanical  seals.  Each 
one  of  these  arrangements  may  result  in  a  different  particular  form  of  eq.  (1) 
making  it  Impossible  to  derive  a  general  kinematic  formulation  of  the  problem. 
This  shortcoming  can,  however,  be  overcome  by  noting  that  the  nutation  angle  y 
in  any  practical  seal  is  very  small.  Hence,  for  small  y  the  spin  ^  in  a  general 
joint  as  shown  in  Fig.  3  can  be  expanded  in  the  form 


(}|  =  ,()q  +  (jilY  +  +  ...  (2) 

where  (|i^  =  are  general  periodic  functions  of  B,  and  B  =  -'ll.  For  y  =  0  '"'e 

have  i|)  =  B,  hence 


=  B 


t  i 


Differentiating  eq.  (2)  with  respect  to  time  yields 


•  •  ^01  *  *  V0O  *  * 

^  =  B  +  —  By  +  (tiY  +  -  6y^  +  2(ti  yy 

dp  dp 


i  I 

*r  ^  - 


The  transmission  law  of  the  joint  is  thus  given  by 


i  3'f'i  I  o  ; 


i.  A 


For  Y  =  0  any  universal  joint  results  in  T  =  1.  Hence, 


L  I 


4,,  =  0 


and  the  transmission  law  can  be  written  as 


3*2  2  yy 
T  =  1  *  — Y  2.1.2^+  ... 
3B  $ 


For  small  perturbation,  the  order  of  y  is  the  same  as  that  of  y»  and  eq.  (5) 


a 


takes  the  form 

T  =  i  =  1  +  0(y'3 
6 


Hence,  for  any  practical  mechanical  seal  where  y^  «  1  equation  (6)  gives 
the  transmission,  T,  accurately  enough,  by  T  =  1. 


*•*  N*  V 

tn 


V.v 


In  a  Constant  Velocity  Joint  the  result  T  =  1  is  accurate  independent  of  y 
This  special  case  is  characterized  by  the  lack  of  preference  to  the  order  and 
direction  of  the  two  perpendicular  tilts  of  the  inner  ring,  and  is  typical, 
therefore,  to  such  cases  where  the  pins  of  Fig.  2  are  either  omitted  or  are  not 
in  effect.  The  first  case  corresponds  to  a  flexible  support  that  consists  of  a 
metal  bellows.  The  second  case  occurs  when  the  friction  in  the  elastomeric 
secondary  seal  is  sufficient  to  prevent  rotation  of  the  flexibly  supported 
element. 

As  another  example  let  us  examine  the  case  of  a  Hooke  Joint.  The  equation 
of  kinematic  constraint  for  this  particular  joint  is  (see  p.  272  in  ref.  (7)). 

tan  (J)  =  tan  B  cos  y  (7) 

Differentiating  with  respect  to  time  and  substituting  tan(^  from  eg.  (7),  gives 
after  some  algebra 


«  • 

,  Bcosy  -  YSinYsinBcosB 
^  =  1  -  sin'Bsin^Y 


For  small  nutation  angles,  «  1,  we  may  use 


sin  y  =  y 


so  that 


cos  Y  =  1  -  ^ 

(1  -  sin*Bsin*Y)~‘  =  1  +  Y*sln*B 


(8) 


Substituting  these  relations  in  eq.  (8)  and  neglecting  terms  of  order  higher 
than  Y^,  will  finally  give 


(9 


T  =  T  =  1  +  (sin^8  -  i  )y^  -  sinB  cosB 
6  B 

Comparing  corresponding  terms  in  eqs.  (9)  and  (4),  we  have 


and 


d^2 

ii” 


=  sin^B  -  i. 

2 


which  agrees  with  the  fifth  term  in  (4)  that  is 


it  =  -  i  sinBcosB 
2  2 

The  particular  form  of  eq.  (2)  for  a  Hooke  Joint  is  therefore 


B  —  sinB  cosB 
2 


(II 


Recalling  that  B  =  -Ip  and  hence  B  =  -<P  we  may  conclude  that  for  any  praC' 
tical  mechanical  face  seal,  where  y^«l,  eq.  (6)  gives  the  spin  accura¬ 
tely  enough,  by 


(1. 
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THE  DYNAMIC  MOMENTS 


The  term  "dynamic  moments"  is  used  to  describe  the  contribution  of  the 
inertia  of  a  body  to  its  behavior  in  the  angular  degrees  of  freedom.  The  rota¬ 
tional  equations  of  motion  of  a  body  are  formed  by  equating  the  dynamic  moments 
with  the  "applied  moments"  that  are  contributed  by  external  forces  acting  on  the 
body.  The  correct  formulation  of  the  dynamic  moments  that  act  upon  the  flexibly 
mounted  element  of  the  seal  is,  therefore,  essential  for  the  two  basic  arrange¬ 
ments,  the  kinematics  of  which  was  analyzed  in  the  previous  section. 


The  general  form  of  the  dynamic  moment  vector  of  a  rigid  body  expressed  in 
a  moving  reference  can  be  found  in  several  texts,  e.g.  [?]  and  is  given  by 


^  ^  ’'go''  "^o 


(12) 


where  L  is  the  relative  angular  momentum  vector  of  the  rigid  body  defined  as 


{L}  =  fll  M 


(13) 


and  X  is  the  absolute  angular  velocity  vector  of  the  body.  The  vector  is 
the  rotational  velocity  of  the  reference  system  accelerating  at  a^,  and  is 
the  location  of  the  center  of  mass  of  the  body  (in  our  model  r  =0).  The 
absolute  angular  velocity  X  of  the  body  is  given  by 

t  =  u  +  (14) 

where  SI  is  the  angular  velocity  vector  of  the  body  relative  to  the  rotating 


reference. 


In  the  seal  model  shown  in  Fig.  2  the  body  is  the  inner  ring  and  the 
rotating  reference  is  the  xyz  reference.  Hence,  due  to  the  kinematic  constraint 
ft  is  always  along  the  z  axis  and  by  definition  is 

=  i^z  (15) 

•  • 

The  spin  ({i  is  related  to  the  precession  ij)  by  the  equation  of  the  kinematic 
constraint,  which  for  small  nutation  y  is  given  in  eg.  (11). 

Flexibly  Mounted  Stator;  The  angular  velocity  of  the  rotating  reference  for 
this  case  is  (see  Fig.  2) 

A  •  ^  A 

=  YgX  +  siny^y  +  cosy^z  (16) 

Where  the  subscript  s  is  used  to  indicate  the  stator  as  the  flexibly 
mounted  element.  Substituting  eqs.  (16),  (15)  and  (11)  in  eq.  (14)  we  have  for 
the  angular  absolute  velocity  of  the  stator 

A  •  A 

siny  y  +  ijl  (cosy  -  l)z  (17) 

O  O  O  ^  O  d 

Hence,  by  eq.  (13)  the  relative  angular  momentum  vector  L  is 

^s  ""  ^^s’^  I'l's^iny^y  +  I^'l’gCcoSYg  -  l)z  (18) 

where  I,  is  the  polar  moment  of  inertia,  and  I  =  I  =  I  is  the  transverse 
z  X  y 

moment  of  inertia  of  the  flexibly  mounted  element. 
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Using  eqs.  (18)  and  (16)  in  eq.  (12),  recalling  that  =  0, 

and  that  we 

are  dealing  with  small  angles  y  so  that  cosy  =  1  and  siny  =  y 

s  s  s 

are  valid 

i 

approximations,  we  have  the  dynamic  moments  in  the  form 

..  ,2 

T  =  I(y  -  Ip  y  ) 

X  '■ 's  ^s's 

(19a) 

4 

' 

••  •  •  .  1 

(19b) 

(19c) 

.4 

As  can  be  seen  from  eq.  (19c)  the  dynamic  moment  is  of  order  y*  and  hence, 
can  be  neglected  in  any  practical  seal. 


Flexibly  Mounted  Rotor i  In  this  case  the  outer  ring  in  Fig.  2  which  represents 
the  shaft  has  an  angular  velocity  u.  This  velocity  when  added  to  the  relative 
precession  Ip  of  the  rotating  reference  xyz  gives  the  absolute  precession  of  the 
rotor,  ip^,  in  the  form 

Ipj.  =  Ip  +  u)  (20) 


where  the  subscript  r  is  used  to  indicate  the  rotor  as  the  flexibly  mounted 
element. 


From  eqs.  (11)  and  (20)  we  have 

•  • 

Ip  =  u>  -  Ipj.  (21) 

The  angular  velocity  o)^  of  the  reference  XYZ  is  given  by  eq.  (16)  where  the 
subscript  s  is  replaced  everywhere  by  the  subscript  r.  Similarly,  the  absolute 
angular  velocity  of  the  inner  ring  is  given  by  eq.  (14)  where  0  is  given  in 
(15)  and  Ip  in  (21). 


»  V  -  V.v  vjui.y  *.*"-V-*  -*■'  'J*.*^.'  ^  *  J  *M  »v:'  «■*  ' 


Hence, 

^  A  .  A  *  A 

=  Yj,x  +  ‘l^j.sinYj.y  +  [ll)j.(cosYp  -  1)  +  uJ  z 

(22) 

-> 

The  relative  angular  momentum  vector  L  is 

Lr  =  lYpX  +  I'l^j.sinYj.y  +  I^t'l'j.CcoSYj.  -  D  +  ailz 

(23) 

and  the  dynamic  moments  have,  by  (12),  the  form 

•  •  •  • 

T  =  I(y  ~  ull)  Y 

X  ^ 'r  ^r’r'^  z  ^r’r 

(24a) 

Ty  =  KHi^Yj-  2liij,Yi.)  -  I^WYr 

(24b) 

Tz  =  -IzCi^Vr  ^ 

(24c) 

Here  again  T^  is  of  order  and  can  be  neglected  in  practical  seals. 
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SUhtWRY  AND  CONCLUSION 


The  kinematic  model  of  mechanical  face  seals  was  presented.  Two  basic  seal 
arrangements  were  considered.  These  are  the  flexibly  mounted  stator  and  the 
flexibly  mounted  rotor.  The  kinematic  constraint  provided  by  the  antiroation 
locks  in  the  first  arrangement  or  by  the  positive  drive  devices  in  the  second 
was  shown  to  be  similar  to  that  of  a  universal  joint.  It  was  shown  that  inspite 
of  the  numerous  variations  of  antirotation  locks  and  positive  drive  mechanisms 
found  in  mechanical  seals,  it  is  possible  to  present  the  equation  of  kinematic 
constraint  in  the  simple  form 

j,  =  -III 

This  unified  relation  is  the  result  of  the  very  small  nutation,  y,  in  prac¬ 
tical  seals,  and  is  accurate  to  an  order  y*  »^ere  y  «  1. 

The  simple  general  form  of  the  equation  of  kinematic  constraint  enables  to 
derive  the  dynamic  moments  that  act  on  the  flexibly  mounted  seal  element.  These 
moments  are  presented  in  eqs.  (19)  for  the  case  of  a  flexibly  mounted  stator, 
and  in  eqs.  (24)  for  the  case  of  a  flexibly  mounted  rotor.  In  both  cases  the 
dynamic  moment  T2  which  is  the  axial  component  of  the  moment  vector  was  found 
negligible.  The  two  other  components,  namely  Tx  and  Ty  depend  on  the  transverse 
moment  of  inertia,  I,  in  the  the  case  of  a  flexibly  mounted  stator,  and  on  both 
the  transverse  and  polar  moments  of  inertia,  I  and  Iz  in  the  case  of  the 
flexibly  mounted  rotor.  The  contribution  of  the  polar  moment  of  inertia  in 
seals  with  flexibly  mounted  rotor  alters  the  dynamic  moments  Tx  and  Ty  as  com- 


pared  to  the  flexibly  mounted  stator  case.  This  is  equivalent  to  altering  the 
inertia  of  the  flexibly  mounted  element  and  may  affect  the  dynamic  behavior. 


The  analysis  presented  in  this  paper  assumed  no  more  than  two  pins  as  a 
representation  of  the  constraint  provided  by  the  antirotation  locks  or  positive 
drive  devices.  If  three  or  more  units  are  effective,  then  the  inner  ring  is 
actually  “locked**  and  is  unable  to  track  angular  misalignment  of  the  rigidly 
mounted  element.  Such  a  condition  can  be  avoided  by  limiting  the  number  of 
antirotation  locks  and  positive  drive  devices  in  a  seal  to  two  units  at  most. 
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A/  =  dimensionless  moment.  M*fSr,] 

M\,  =  dimensionless  moment  due  to  stator  initial  misaliKiiment 
A/y,  =  dimensumless  moment  due  to  rotor  runout 
m*  =  stator  mass 
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I*  dimeiisioidrsH  piesMiie. ///.S 
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the  lelative  misalignment  betireen  the  mating  faces  is  ^iven  as  a 
fninlion  of  rotor  innout.  assembly  misalifenment,  design  parame¬ 
ters,  and  opeiation  fonditions.  An  expression  is  provided  for  the 
erifnal  rotor  lunout  above  which  the  seal  will  fail  due  to  face 
contact.  .Althoufrh  the  analysis  is  based  on  small  perturbation  as¬ 
sumption,  it  is  shown  to  be  valid  in  many  practical  cases. 

INTRODUCTION 

t'vKlentes  of  dyiianiu  problems  in  mechanical  face  seals 
exist  /dr  more  than  two  cfecades  {/).  Vet,  most  of  the  lit¬ 
erature  on  fluid  sealing  concerns  steady^siale  conditions. 


3*  =  fa<  e  tonittg 

3  -  tiimeusionless  toning. 

7*  =  relative  misalignment 

7  -  tlimcnsionless  misalignment.  7*r, /(.'„• 

7..  =  tfiineiisionless  relative  misalignment  caused  l)N  rotor  iiin- 

out  alone 

7,  =  dimensionless  rotor  runout 

=  dimensionless  criiieal  rotor  runout 
7.  =  dimensionless  lilt  of  the  stator 

7«,  =  dimensionless  stator  initial  misalignment 

7a  =  dimensionless  steady-state  stator  response  due  to  7.,  alone 

7.,  =  dirnensiofiless  steady-state  stator  response  due  to  7,  alone 

6  =  angular  c(M)rdinate 

p  =  viscosity 

%  =  damping  ratio 

c|f  =  precession 

(u  =  shaft  angular  velocity 

Subscripts 
rr  =  critical 

d  -  dynamic  (vehn  iiy)  ef  fect 

e  =  elastomer 

/  Ihiid  him 

i  iiiiiei  radius 

M  =  mid  radius 

o  =  outer  radius 

r  =  nMor 

s  -  siai<»r.  or  flexible'  siip|K)ri 

\p  =  springs 
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()iil\  lecrnllv  lias  dviiatiiits  Ix-nmif  iln-  siiliiiii  <>l  itii-<ii«‘i 
u  .il  iiiM  slin.iliims  .111(1  tnoic  i  inoious  cxim'I  iiiiciil.ilioii.  Rcl- 
fifinv  {!)  (dn-is  llif  lilfialiiif  iiDlil  I’.WI  wiiilc  ReK.  (2)- 
(  /)  presenl  more  lereiil  results  on  the  sul))e(  t. 

The  t).isi(  se.il  model,  .is  shown  in  1.  has  one  hxed 
( ornponent  th.it  is  ri^idh  mounted  either  to  the  shaft  or  to 
the  housing  (m  Fi^^.  I  the  ri^idh  mounted  (omponent  is  the 
rotor).  File  other  (  iimjxmenl  is  tlexihiv  inotinted  to  allow 
tr.u  kiti^;  ol  \  .n  iotis  ims.diKtmtents.  I  hese  tmsali^nttientsr.m 
Ik'  the  lesiilt  ol  m.imil.ii  tut  ing  ot  asseitthh  toletaiues.  I  hits. 
Ill  the  seal  ol  Fit;.  1,  lot  ex.imple.  the  rolot  m.n  have  .i 
(Ctt.iin  .ixi.il  Minoiit  ,md  the  st.ttoi  too.  m.iv  Ih-  misaligned 
III  llie  sli.ill  due  lo  impel  lei  lions  m  Ms  Ilexihit'  siip|Hiil.  In 
l.ul.  It  IS  \eiv  iinlikelv  lo  li.oe  .iiiv  pel  le<  I  .ilignmeni  .it  .ill 
III  pi.Klu.il  se.ils  In  opei.ilion.  the  loioi  wolihles  liei.iiise 
ol  Its  I  mioiil  iiid  this  motion  is  1 1 .iiismilled  vi.i  the  Ihiid 
film  lo  die  St. 1(01 .  I  he  dvii.imii  ( li.ir.ii  lei  isIk  s  ol  die  svsiem 
iiumielv.  stillness  and  d. imping  of  liolli  the  fluid  him  and 
Ilexihle  support,  .is  well  .is  the  ineili.i  ol  the  si. nor)  will 
del  ei  mine  the  dvii.imii  In'liavioi  ol  the  se.il,  Fhe  stator  <  .m 
hei  ome  (  omplelelv  iinstahle  or  il  t  an  it  .n  k  the  rotor  runout 
in  .1  d\ n.iinii  alls  st.ihle  mode.  Howevei ,  even  in  this  stable 
mode,  the  i oinhin.iiion  ol  traeking  amplitude  and  phase 
slnli  1 ,111  he  su(  h  ih.il  the  relative  (Misition  between  stator 
and  Kitor  results  in  seal  lailtire.  Siu  h  a  lailtire  lati  he  either 
due  to  excessive  leakage  caused  hv  huge  gaps,  or  due  to 
loi.il  l.ue  (ont.iit  in  high-speed  nomotitacting  seals. 

Iti  ptevious  woi  ks  on  seal  dviiamics.  e.g.  (5).  (6),  some 
simplilviiig  assumptions  were  made  which  restricted  the 
results  hut  neveriheless  g.ive  good  insight  oh  the  prohletti. 
File  I  esli  icliiig  .issutnplions  iiu  hide  neglei  t  oF  various  lltiid 
him  el  lei  IS.  neglei  I  ol  rotor  and  slalot  niis.ilignnients,  atid 
iiegled  or  inioinplele  iiilormalioti  on  the  dviiamii  piop- 
ei lies  ol  die  Ilexihle  siippoi  1  In  this  paper,  a  dvnaiiiii  aii.il- 
V  SIS  ol  ( oned-I.K  e  se.ils  is  presented  wliii  h  lakes  into  avi  iiiinl 
’ll. IIIV  lailors  di.ii  .died  die  seal  dvn.uniis.  1  hese  ladors 
.ire  the  roioi  i  iinoiil.  initial  tiiisalignnieni  ol  the  stator  l>e- 
loie  .iit.K  hnieni  to  the  rotor,  coning  ol  the  slalot,  slilTness 
.md  d.iniping  of  the  Ihiid  film,  and  siif  lness  and  damping 
o!  the  Ilexihle  support.  I  wo  assiiinplions  will  1k'  made,  how¬ 
evei.  to  enable  ilosed-lorm  analvliial  solution,  ft  will  Ik' 
.issniiied  dial  iDihe  seal  is  onlv  slightly  pei  liii  hed  Iromlhe 
position  ol  p.irallel  f.iies.  and  (‘J)  that  the  sl.ilot  anipliliide 
IS  small  enough  lo  prevent  slippage  of  the  secondarv  seal. 


SPRING  PRIMARY  SEAL  RING 


I  he  lust  ,issumpiion  seems  ,il  Inst  lo  Ik-  veiv  lesiiidive 
llowevei.  It  will  he  shown  hv  ,i  niiinerii.il  ex.iniple  that  this 
IS  die  lase  in  most  pradii.il  applications,  1  he  seioiid  as¬ 
sumption  is  valid  lot  manv  practical  cases  and  permits  rep¬ 
resentation  of  elastoniet  i(  secondarv  seals  hv  relalivelv  sim¬ 
ple  lo  handle  dvnamic  coel  licients. 

.'Mthoiigh  the  present  analvsis  is  performed  lor  a  coned- 
laie  seal  with  an  elastomeric  O-ring  secondaiv  seal,  the 
results  will  he  presented  in  lerms  of  general  dynamic  coel- 
li(ienls  and,  hence,  are  applicable  lo  other  nonconlacling 
se.il  ivpes  as  well, 

SYSTEMS  OF  COORDINATES  AND 
GENERALIZED  FORCES 

I  he  seal  model  and  v  at  ions  i  oordinale  systems  ate  shown 
111  Fig.  2.  Fhe  inertial  leleience  XYZ  is  fixed  in  spat  e  and 
its  axis  /  coincides  with  the  shall  axis  of  rotation.  Fhe  system 
xvt  delinc-s  the  relative  position  of  the  stator  with  respect  to 
the  inertial  releretue.  I  bis  svsiem  is  located  in  the  stator 
plane  hut  is  Free  lo  rotate  so  that  axis  v  always  points  lo  the 
|M)ini  ol  maximum  dislani  c'  Ironi  the  plane  X>’,  and  axis  x 
is  always  parallel  to  the  ,V)'  plane.  Fhe  angle  ili.  between 
axes  X  and  x,  is  the  precession  and  the  angle  y*  is  the 
niilalion  defining  the  .ingular  tilt  of  the  stator  with  respect 
to  the  inertial  reference.  Fhe  system  x,y,i,  is  fixed  to  the 
rotor  and  rotates  with  it  at  the  angular  velocity  to  about  the 
7.  axis.  Axis  x,  is  always  parallel  to  the  XY  plane  making  an 
angle  col  with  the  X  axis.  .Xxis  y,  always  points  to  the  point 
of  minimum  distance  from  the  XK  plane.  The  angle  yf 
represents  the  runout  of  the  rotor.  Because  the  separation 


2 


Fig.  1— SctMfiMtle  of  a  radial  laea  aaal 


Fig.  2— Saol  modal  and  eoordinala  ayalama 


Slai)ilil\  rhr«’siii>l(l  and  Sit‘ad\ -Stale  Response  «>1  Noneoniailin^  (ained-raie  Seals 


(.  tielween  die  sialnr  and  roioi  Is  \ei  v  small  (a  lew  ink  I'o- 
inetei  s  in  prat  lit  al  applieations)  -y*  anti  -/?■  are  alsti  very  small 
anti  t  an  lie  irealetl  as  vet  ltirs.  Hent  e,  die  relative  misalign- 
menl  lielween  die  slaliir  anti  rotor  is 


Hgiire  S  shows  die  relalive  posiiitm  iK-lween  die  stator 
(seal  ring)  and  rolor  (seal  seat).  I  liis  relalive  position  is 
ilestril»'il  liv  .1  new  t  tioi  tliiiale  system  12;).  I  liis  system  is 
tree  to  rol.ile  wiilim  the  sialoi  plane  so  that  axis  1  (ahoui 
whit  li  the  slaloi  lelalive  lilt  -y*  lakes  plate)  is  alw.iys  paiallel 
to  the  loloi  pl.iiie  anti  axis  2  always  points  to  the  ptiint  of 
m.iximnm  him  ihit  kness. 

I  he  y.niotts  loorilinale  systems  anil  a  yeilor  represen¬ 
tation  til  the  tills  -y,  -y,.  anil  "y,  alioni  the  axes  I,  .y.  anil  x,. 
I espei  lively ,  ate  shown  in  Fig.  1.  I  he  .nigle  tji  is  measmetl 
lioni  the  pl.ine  /)  to  the  plane  ( ),(  t,.VI{  in  Fig.  2.  The  angle 
It  IS  measiii  eil  ini  he  sl,iloi  pi, me  1 1  om  axis  2  to  I  he  line  OJt. 
Iletiie.  the  loi.il  film  ihit  kiiess  me.isineil  helyveen  points /t 
.mil  H  in  Fig.  2  is 

h  =  C  -1-  ■y*rfose  t  (l*(r-i,)  12] 


e  =  4.  -  tb,  [:M 

anti  tbi  is  the  angle  hetween  .ixes  \  anil  1  or  alternatively 
between  axes  .mil  2.  Fi|iialion  (2|  lakes  attiiitnl  ol  stator 
I'omng  as  shoyvn  in  Fig.  If  yvhere  3*  is  the  ttining  angle. 

I  he  wohhiing  moiitm  ol  the  rolor  due  to  its  riinotit,  anti 
the  lelalive  rnis.ilignnieni.  y.  between  siatoi  and  rotor  pro- 
tliK  e  a  system  of  fiiiid  him  loites  and  moments  which  is 


.ipplietl  to  the  stator.  I  he  flexible  sup])t)rt  ol  the  siaioi 
which  consists  of  spi  ings  and  a  secondary  seal  (tisiialh  an 
elasiomerit  ring)  proihii  e  .mothei  system  ol  Ion  t-s  .mil  mo¬ 
ments  that  an  tm  the  si.iioi.  In  lineal  systems,  the  sum  of 
the  applied  forces  in  ilegree  ol  freedom  j  title  to  a  ilisiiirh- 
ance  in  tiegree  of  li  eetlom  /  can  be  expressetl  in  I  lie  general 
lol  III 

2^/  =  I'll 

where  F,  is  a  geneiali/eil  force  (a  force  or  a  moment)  and. 
K,,  anil  I), I  are  stillness  ,mil  ilamping  coefficients,  respet- 
livelv.  In  our  case  (where  lor  small  perturbations  the  system 
can  be  consklcretl  linear),  both  the  Huid  him  and  the  flexible 
support  itintribnle  stillness  and  damping.  Hente,  the  siilT- 
ness  K  IS 

K  -  Kf  -1-  K,  (.">a| 

where  the  support  siiliness  K,  includes  the  contributions  of 
the  springs  anti  the  elastomeric  secondary  seal.  Thus. 


T  he  damping  1)  is 


K,  =  K,,,  +  K, 


1)  ^  D,  +  D, 


where  in  many  practical  cases  the  only  source  for  the  sup- 
(Mirt  damping  D,  is  the  elastomer,  in  which  case  /),  =  D,. 

Fluid  him  stiffness  and  damping  coefhcients  for  a  coned- 
f'ace  seal  were  analyzed  in  (7)  and  are  summarized  in  Ta¬ 
ble  I.  These  coefhcients  are  conveniently  calculated  in  the 
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SEAL  SEAT  -< 


Fig.  3— Ralativ*  position  bstsissn  tool  ring  and  soal  saal 


Fig.  4— Vaolor  rapcaaantatton  of  tho  Hlta  y,  y^  yt 
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T  \BI  f  I  —  XOMIIMI.VSKISAI  Kl  (111  FiI.M  SlIffNISS  .AM>  OAMflM.  (  .<  fH(  IKM  S  K/  AND  D/. 


1  =  1 

/  =  ^ 

/  =  3 

0 

1 

0 

1  -  L' 

1 

0 

0 

1) 

/  -  :i 

K 

0 

(» 

1) 

A-nHj;,, 

(  I  -  /{,)  l{„ 

•J  +  P(l  -H.) 

Oill  +3(1  -«,))  -  2  '  ~-- 

*  '  2  +  3(l-K,) 

3  '(l 


1  -R, 
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loonliii.iu-  I2!l  of  Fi^.  A  and,  lor  the  small  pertur- 

haiioii  t.isf,  .irf  liiuMi  It  should  be  noted,  however,  that 
the  thnaniit  loelluients  presented  in  Table  1  were  found 
s.ilid  eveti  bevond  the  lanKe  of  stnall  pertuibations  (7). 

()\  ii.iinu  piopei  IK  s  o(  e(asioineri(  ( bring  seoindai  v seals 
wcic  iioeslig.iled  ;ii  (.S'),  In  the  .ibsente  ol  slippage,  their 
sill  I  ness  ,tnd  d.unping  <  oel  Ik  ienis  <  an  be  ex|)iessed  in  the 
lot  in 


l\.  ~  Aio" 

.nid 

f  1  /> 

I),  -  au) 


lb) 

Ib.-O) 

Mui'Km  ?bi 
-  hMtf  in  ipac* 


ulieie  the  oielliiienls  .i  H.a.h  ilejH'nd  on  lailors  sueh  as 
iii.ileri.il  loniposiiion,  lein| 'eralure,  ring  dimension,  pre- 
lii.id  eu 

In  ohiaining  the  system  of  fortes  and  monienis  .ipplied 
bv  the  Ilexible  support,  it  is  useful  to  observe  the  various 
stages  of  seal  assembling  as  shown  in  Fig.  .a.  Belore  any 
■itiat  hineni  t.ikes  plat  e  [Fig.  .a(a)|,  the  roltir  has  its  runout 
y,  and  the  si.iior  has  a  teriain  initial  misalignment  y„-  Uur- 
iiig  the  .isseniblv.  (he  st.iior  is  pressed  against  the  rtilor 
I  Ins  I  .III  Ih’  dest  I  iIh’iI  bv  In  si  tilting  the  stator  bv  the  .imoiinl 
y.,  |Fig  a(bl|  whiili  ie<|uires  ,i  nionieni  .\/v,  given  bv 


(cl 

“r,-K,H  ’'r 
rotating  momont 


Fig.  5 — Ooocriptlon  ot  Initial  mloallgninanta  and  attachmant  aaquanea  of 
Iho  saal  comptMianta. 


’  Iv.iiy, 


M'l 


ttlieie  /v.ii  is  the  .ingular  stillness  ol  the  Ilexible  sup|>orl 
(see  .\ppendix  I).  I  be  inili.il  lilt  of  the  st.itoi  was  seleeied 
.irbilrarilv  in  the  op)H>site  diieilion  ol  the  inertial  axis  X 
without  losing  the  generaliiv  of  the  .inaivsis.  Sinte  y.,  is 
fixed  ImiiIi  in  magniliide  and  diretiion.  the  moment  Af\,  of 
Fit).  |(>|  IS  also  (ixeil  in  magnitude  and  diiet  tion  (whit  h  is 
along  ihe  .V  axis). 


Final  allathnieni  of  the  seal  mating  faces  in  the  assem¬ 
bling  pnness  |  Fig.  ")(( )|  results  in  an  additional  moment  that 
is  re(|nired  lo  till  the  stator  by  the  amount  of  the  roior 
rnnoiil.  Hence. 

iW,.  =  K,ny.  [7] 
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Siahilitv  I  lirt-shoUi  and  Sieady-Siate  Response  of  Noncotuactiii)'  Cx>ncd-Face  Seals 


The  moment  A/.„  like  the  runout  y,.  has  a  fixetl  magnitude 
hut  is  rotating  at  the  shaft  speetl  o).  and  is  tlirected  along 
the  axis  x,  (Fig  -1). 

I  he  moments  acting  on  the  stator  during  operation  atid 
resulting  from  the  flexible  support  are  given  in  the  system 
xyz  of  Fig,  1  in  the  form 


=  -A.ii7,  -  /Tii7, 

[«al 

.V/„  ^  -//,,, iji7. 

|8b] 

Fluid  film  tnoinents  expressed  in  terms  of  the  fluid  film 
stiffness  and  damping  coefficients  (7)  are  given  in  the  ref- 
ei  elit  e  svstem  f'J.'l  liv 


Mi\  -  -A'/ii7  -  Di  i[y  (‘fa| 


I  he  complete  svstem  of  moments  acting  on  the  stator  is 
shown  in  the  vec  tors  diagram  of  Fig.  (i.  With  tlie  aid  of  this 
figure,  there  is  no  cliftic  lilts  in  clesc  rihing  tlie  tnoments  in 
.iny  t  onvenient  reference  svstem  like  the  inertial  svstem  XiV. 
or  the  rotating  svstem  xyi. 

1  he  net  axial  force  along  the  /  axis  is  gi\en  hv 

=  -  A.  ,,/  -  l)„/.  [1(1] 

SEAL  RING  STABILITY 

In  the  following,  cliinetisional  paratneters  will  lie  indi- 
lalftl  hv  ati  asterisk  to  distitiguish  them  f  rom  their  corre¬ 
sponding  dimensionless  sallies.  .\lso,  the  suhscript  II  will 
Ik-  omitted  frciiii  the  angular  stillness  and  damping  coef¬ 
ficients.  ffence,  whenever  A  and  1)  will  Ik'  used  without  any 
index  tfies  will  correspond  to  angular  degree  of  freedom. 
Onis  the  subsc  i  ipi  will  be  retained  for  the  axial  degree 
of  fic-edom. 


^  ,  I 

/  ^  Ml, 


:  V'  ' 

\. 


r  T'  Vs 

I  ^ 


-  ^ 


^  M 
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'm.. 


Pig.  Vaelor  r»pfM«n>rtloo  of  IHta  and  momanto 
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T  he  equation  of  motion  in  the  axial  degree  of  f  reedom 
is  by  [10]  and  [68]  of  Appendix  2: 


ml  +  D  m?.  -I-  KmI  =  0  [II] 

This  equation  is  linear  and  uncoupled  with  the  equations 
of  motion  in  the  angular  degrees  of  freedom.  Its  solution 
is 


/  =  zy 

resulting  in  the  characteristic  equation 

mk’  +  +  A:,:,  =  0  [12] 

The  condition  lor  stable  operation  in  the  axial  mode  re¬ 
quires  that  all  the  coefficients  of  [12]  will  lie  positive.  Fhe 
dimensionless  mass  m  and  support  damping  /An.  are  alwavs 
[jositive.  The  damping  /f  ic  is 

On  =  /),n  +  D/n 

and  since  lif  w  >  0  (see  I  able  1 ).  /Ais  is  also  alwavs  |H>siti\e. 
Hence,  in  order  to  satisf  y  the  stability  tondition.  it  remains 
to  iec|uire  that  An  >  0.  Substituting  A/ss  from  Fable  1  in 
[.'ra],  we  have 


W(0,.  -  Ai;  +  A,n  >  0  (13) 

C-ondi  I  [13]  is  satisfied  whenever  p  =  0  (corresponding 
to  a  face  seal).  With  coned  face,  axial  stability  is  always 
obta.  ed  when  the  film  converges  in  the  direction  of  pres¬ 
sure  drop. 

T  o  analyze  the  angtilai  stability,  the  linear  equations  [66] 
and  [67)  of  Appendix  2  are  used.  For  the  transient  lesfKinse, 
only  the  homogeneous  part  of  these  equations  has  to  be 
considered,  namely 


h\  +  f>‘Vv  +  Kyx  +  ^0/lrr  =  «  [ITa] 

h)  +  Oin  +  Ayi  +  -D/yx  =  0  [14b] 

Using  complex  notation,  the  relation  between  the  stator  tilt 
7,  and  its  components  yx  and  7,  can  be  represented  in  the 
form  (see  Fig.  4) 


7.  =  7x  +  nr 


[1.5] 


witere 


1  - 

Miiltipiving  F.q.  [  I4h]  by  1  and  adding  it  to  Eq.  [14a],  noting 
that 


n~  =  »7v  -  7r 


() 


I.  (iRftS  \\l>  I.  K  IMo\ 


V  l('l«  Is 


Mtr  [Kir  m-^)  >  0  jjh  1 


h  +  ih-  +  Ky~  -  i7,ihy'  =  " 

I  lie  loniplex  leini  -(-/>/7,  in  K<|.  [lt)|,  indiued  by  hv- 

diiidsM.niui  el  lei  ts  in  die  tiuid  film,  represents  a  transverse 
iiioineni  uliitli  le.ids  the  displatement  y,  b\  it''.*  radians, 

I  he  MiluiKin  III  I  (|.  I  l(i]  is  III  the  liiini 

y.  -  y. IKia) 

and  ihei  omples  ■  li,ii  ,ii  tei  islii  e<|natiiin  is 

IK-'  r  l)K  ^  K  -  iH  -  0  |I7| 

where 

H  -  n,/'!  . 

MiiliipK in^i  K().  1171  hs  its  <<injuj»a(e  lorni  gives  the  ehar- 
aiieristu  eipiaimn  Im  the  angular  mode  in  its  real  iorm 

/■'a'  ^  21  Dk''  +  VJIK  +  0-’)A- 

1  InJ 

+  2I>KK  +  +  «“)  =  0 


411)^  (A-  +  li-hKl)-  -  lli^)  >  0  (L’ld) 

.Since  /  and  1)  are  positive  and  lulhllment  ol  condition  120) 
means  that  A  is  positive  too,  it  remains  to  require  that 

KD^  -  IB~  >  0 

Hence,  the  angul.ii  siahilliv  threshold  lorresponds  to  the 
<  ase 

KD-  -  Hi-  -  0  (22) 

.Siihstituting  H  =  /)^/2  and  learranging  Eq.  |22]  gives  the 
dimensionless  c  ritic  al  moment  of  inertia  in  the  form 

/„  =  ^(A;  +  A,)^l  +  [23a] 

In  Ref.  (6),  an  empirical  stability  threshold  was  found  for 
the  case  /X  =  0  from  a  vast  number  of  numerical  solutions 
without  the  restriction  of  small  perturbations.  The  critical 
moment  of  inertia  in  (6)  has  the  form 

h,  =  4(A,  +  RfA/)  [23b] 


I  he  neiessarv  cotidiliou  for  stability  requires  that  all  the 
iiicihi  leiits  m  I  IS|  will  he  positive.  I  he  ditiieiisionless  mo- 
iiiciii  ol  iiu  ili.i  /  IS  .ilw.iNs  positive,  l•'lom  [."ihl  and  fable  1, 

II  i.iii  he  si'cn  lhai  l>  ■  II.  Heme,  angular  stahilily  re<|uires 
ill. II 


2/A  +  1)^  >  (1 


[19al 


.iiid 


2/;A  X)  [191)1 

(.oiidilions  [19a]  and  [I9h]  are  salislied  whenever  A'  >  0. 
I  liiis.  trom  Eq,  [.aa]  and  Table  1,  the  necessary  condition 
|oi  .iiigular  stahilils  is 

■nil’,.  -  /’,)(P/?,  -  DA.'f  +  A,  >  0  [20] 

Defining  a  critical  coning  angle  3,,  =  1/W,.  it  is  seen  from 
[20]  that  foi  /’„  >  /’,  .mgular  stability  is  assured  whenever 
3  >  3/>.  In  Kef.  (7).  .Ill  optimum  cotiing  atigle  for  maximum 
angular  stillness  K/  w.is  found  in  the  lortn  p.yv  =  2//f,(l  -H,). 
It  is  clear  ih.il  this  3../.r  is  larger  than  3'>  •<>'<1.  hence,  an 
optimum  coning  angle  assures  the  lullillmeni  ol  condition 
(2U]  whenever  /’„  ■ 

file  Koiith-Hui  w  it/  stability  criterion  for  Eq.  [18]  re- 
t|uiies.  in  addition  to  fullillitig  conditioti  [20],  that  the  four 
following  conditions  are  also  met: 

211)  >0  [2 1  a] 


where  K/  has  the  value  of  the  Huid  film  coefficient  A/n  of 
Table  1.  As  can  lie  seen.  /,,  of  [23b)  is  very  similar  to  that 
of  [23a)  when  //,  =  0,  esjK'cially  lor  narrow  seals  where  R, 
is  c  lose  to  unity. 

Replacing  the  dimensionless  parameters  in  Eq.  [23a]  with 
their  corresponding  dimensional  parameters  (see  nomen¬ 
clature)  yields  a  dimensional  critical  speed  ci)„  in  the  form 


Cl) 


2 

rt 


Kf  -C  A* 


[24] 


Eoi  atigular  stability,  the  operation  sfieed  ii>  should  be  smaller 
than  ci),,.  As  ran  l)e  seen  from  [24],  an  increase  of  the  stiff¬ 
ness  and  reduction  of  the  inertia  increases  the  critical  sjieed 
as  expected.  However,  increasing  the  fluid  film  damping 
Of  reduces  co,,.  This  is  due  to  the  hydrodynamic  effect  (rep¬ 
resented  by  D*)  which  results  in  a  transverse  moment  lead¬ 
ing  the  stator  tilt  by  •ir/2  (see  Eq.  [16]).  Since  K*  is  directly 
pro|)ortional  to  the  pressure  differential  p„  -  p„  and  D*  is 
directly  pro|K)rtional  to  the  viscosity  |i,  it  can  be  seen  from 
[24]  that  high-pressure/low-viscosity  conditions  improve  the 
seal  stability. 

It  should  l>e  tinted  that  A,*  and  D*  in  Elq.  [24]  may  include 
elastomerk  stiffness  and  damping  components  which  are 
frequency  dependent.  In  such  cases,  an  iterative  process  is 
retpiired  in  solving  E^q.  [24]  lor  to,,. 

It  is  of  interest  to  obtain  the  preseccion  iji  of  the  stator  at 
the  stability  threshold.  The  getieral  form  of  the  ex|)onent 
A  in  E]q.  [Ifia]  is  A  =  n  -t-  «[>  but  at  the  stability  threshold 
n  =  0  atid  hence. 
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2I^I)K  +  211)  '  >  0 


(21b] 


SiaMliiv  I  lii«-sli(il(l  .111(1  Sic.idv-Sialc  Kes|)oiis«-  <>l  N<»ru«mla< ting  OoiK'd-Fate  Seals 


Siihsoliiltiig  III  l'<|.  |IS|  and  e(|ii.iliiig  se)).ii .ilcK  le.il  .iiid 
iiii.igiii.il  V  |).iils  III  /erii  gives,  lespei  livelv 

/-'lii'  -  (2/A,  +  /r'jiii'-'  +  (A-  +  //-')  -  0 

.iiid 

/li*'  =  A 

The  steady-state  i  espouse  to  an  initial  stator  misalignment 
7„  alone  is  obtained  by  substituting  y,  -  d  in  Eqs.  (6(i]  and 
|2'i|  l•'V|•  Hence,  using  again  the  complex  form  [l.'t].  we  have 

/r  +  t)y>  +  A'y,  -  i^D/y,  =  A,y„  |S()| 


Irom  which  it  can  Ite  easily  ioiind  that 

li  I  2 

lii.f  “  ~  “  — - 

II  I  +  DJI), 

wlieie  the  index  U  indicates  stability  threshold. 


.Adding  (be  s(|uaies  ol  the  right-  and  lell-h.ind  sides  <il 
K<|.  12h|.  resjtei  lively,  gives  the  transmissihilily  ol  the  svsieni 
in  the  lorni 


STEADY-STATE  RESPONSE 

Due  to  their  linear  n.iture.  the  .iiigul.it  eijii.itions  ol  mo¬ 
tion  imil  .111(1  |ti7|  <an  Ix'  solved  sep.ii.iieb  lor  the  elleils 
ol  the  lolling  luiutioiis  i.iiised  bv  y,  .iiid  y.,,  lespeclively. 
I  he  solution  ol  these  equations  will  Ik-  presented  lor  each 
one  ol  the  two  iorcuig  liuntions  in  die  most  appropriate 
leleieiue  system.  Heiue.  lor  the  lott  ing  Itiiulion  resulting 
lioin  the  rotor  iiinoiit  yet  tor  y,,  yvhidi  is  totaling  at  the 

sh. ili  speed,  the  solution  will  lx-  perlornied  in  the  rotating 
lelereme  v,:.  For  the  loicing  lunition  resulting  Irom  the 
imii.il  iiiisaligntneiit  of  the  st.itoi  -y,,.  which  is  lixed  in  space, 
the  solution  yvill  he  perlornied  in  the  inertiiil  relereni e  XV7.. 

Si.iiiiiig  with  the  elled  ol  rotor  runout  alone  re<|uires 
A/y,  -  II,  .Also,  .It  steady  si.iie  and  in  the  ahseiue  ol  y,,.  the 

si. iioi  piei  esses  .It  aioiislaiit  i.ileiji  Its  tilt  vector  "y,,  caused 
hy  the  lol.iling  yeitor  y,.  has  a  lixed  value  y.  =  -y,,  and  is 
roi.iiuig  .It  the  same  speed  as  y,.  1  he  resultant  relative  mis- 
.ihgiimeiil  veitor  y  (see  Fig.  b)  has,  therelore,  a  constant 
y.itiie  loo.  .iiid  is  lol.uiiig  at  the  same  sjreed  as  y,  and  y,„ 
n.imely ,  .11  the  sli.ill  speed.  Hem e.  ill  -  (Fi  =  Iwhereij/and 
(bi  .ire  diiiieiisionless  paranieters  not  iiiali/ed  by  lo.  I'sing 
F.ijs  |8.ij.  lHb|,  |‘l.i|,  |(lb|.  .ind  |()1 1.  and  substiluting  -y,  =  y„ 

(onstant,  .is  yyell  .is  -y  -  ib  =  0,  we  have 

ly  ~  -  A  y  ,  -  A^yiosp  +  -D/ysiup  +  A.y.cosp' 

[2b] 

(I  -D,y„  -  A/ysinp  -  -D/ycmp  -  A'.yrSinp' 

Fiom  Fig  b.  yve  see  that 

yiosp  '  y„  -  y.iosp'  |.,_l 

ysinp  -  y.siiqi' 

Siibsiiiiiiiiig  |27]  III  (2(i|.  using  F'.qs.  I.bal  and  (.'ib|  and  rear- 
rangiiig.  we  li.iye: 

(A  -  /)y„  =  ^Acosp'  +  /F/sinp' j-y, 

'  ~  ^  [2H] 

=  ^-A'sinp'  ^  Dfcnsp'Jy, 


Since  y„  is  fixed  in  space  and  has  a  constant  magnitude,  the 
res(>onse  y.  =  y,/  yvill  also  lie  fixed  in  space  and  will  have  a 
constant  magnitude.  Thus,  y,  =  y,  =  (land  F.q.  (SO)  yields 


D/ 

taniba  =  - 


where  il»,/  is  the  shift  angle  of  the  vector  y,/  with  resjrect  to 
the  real  axis  A'  in  the  complex  platie  XY. 

In  tiqs.  (29]  and  [31].  the  stiffness  A'  contains  the  Mexihle 
support  stiffness  K,  according  to  [5a].  When  the  support 
has  an  elastomeric  secrtndary  seal,  its  stiffness  K,  depends 
on  the  frequency  of  the  forcing  function.  It  is  important, 
therefore,  to  distinguisii  between  the  A,  values  to  be  used 
in  Eq.  [29]  where  the  fon  ing  functiotr  Irequettcv  is  wt,  and 
in  Eq.  [31]  where  the  fori  ing  function  frequency  is  /.ero. 

From  Eqs.  [29]  and  [31],  the  complete  steadv-state  re¬ 
sponse  of  the  stator  can  now  be  expressed  in  the  form 


Since  y,r  is  lixed  in  spat  e  but  y „  is  rotating  at  the  speed  ui, 
the  resultant  vector  y,  will  also  rotate  at  the  speed  (o  and 
its  magnitude  will  vary  cyi  lically. 

T  he  relative  misalignment  Ix-lween  .stator  and  rotor,  which 
is  the  more  iin|Mntaiit  parameter  for  seal  performance,  is 
obtained  from  1 1]  and  |32|  m  the  form 

y  =  y..  +  y./  [33] 


y.,  =  y„  -  yr 

is  the  relative  misalignment  in  the  case  -y,,  =  0.  T  he  vey  tor 
-y,,  has  a  fixed  value  and  is  rotating  at  the  speed  u>.  Figure  7 


s 


1.  (.KthN  .\M>  1.  I'.ISION 


Fig.  7— Vector  representation  of  tilts  of  the  seal  components 


presents  the  vniHius  till  veetors  in  the  inertial  refercnie 
plane  ,YJ'.  From  this  hgnre,  tvc  .see  that 

y~  ~  yt  +  y~i  +  '2y„y,|^^>s^D^  [Ml 

uherc  t  is  the  time  tne.tsiirefi  from  the  instant  at  which  the 
sector  y,.  passes  over  the  vector  ■y,/. 

From  F,<j.  (iHJ,  it  is  clear  that  the  relative  misalignment  y 
is  also  time  dependent  similar  to  y,.  This  is  because  the 
sedor  y„  is  rotating  wliile  7,/  is  stationary. 

In  order  to  obtain  the  value  of  y„.  F.q.  [27]  is  used  and 
y.iosp'  along  with  7rSuip'  substituted  in  Eq.  [26J  for 
y„  -  (1.  This  gives; 


~/y„  =  ^-A'cosp  +  ^0/sinp^7„ 
D,7„  =  ^-Asinp  -  ^D/cosp^7„ 


l.nes.  lor  example,  loriespoml  to  the  tase  7  -  0  whuh  is 
die  most  ilesiiahle  <ase.  In  pi  at  lice,  however,  the  lelalive 
misalignment  is  always  7  >  0  and  its  exact  value  is  lequired 
in  order  to  enable  accurate  t  alculations  of  leakage,  (rittion 
torque,  etc.  In  the  presence  of  Iroth  an  initial  stator  misalign¬ 
ment,  7„.  and  rotor  runout.  7,.  the  relative  misalignment 
becomes  time  dependent  (see  Ft)  [34]),  which,  in  turn,  re¬ 
sults  in  torque  and  leakage  variations  at  the  shaft  f requeiu  v. 
While  it  is  difficult  to  observe  rapid  leakage  \arialions,  torque 
variations  were,  indeed,  measured  at  shaft  frequencies, 
e.g,  (7). 

A  t<H>-high  relatixe  misalignment  can  cause  seal  lailiiie 
due  to  an  excessive  leakage  or  due  to  IcKal  face  coiil.ici. 
Hence,  it  is  useful  to  investigate  the  relative  misaligmiient 
7  as  alfected  by  the  operation  conditions  and  design  pa¬ 
rameters.  From  F.q.  [34],  it  is  clear  that  in  order  to  lechice 
7,  Ixith  7,/  and  y„  should  be  as  small  as  possible.  From  F.qs. 
|3la)and  (37),  it  isclear  that  ihiscan  Ire  achieved  be  keeping 
both  the  stator  initial  misalignment  7,,  and  the  rotor  ruiioiit 
7,  tninintiim.  Increasing  the  fluid  film  damping.  1)/.  and 
stiffness.  A';,  also  hclirs  in  reducing  Ixtth  7,/  and  y„.  1  he 
dynamic  properties  of  the  flexible  sup|)c)rt  have  a  smaller 
effect  on  7,/  and  y„  hut  reducing  the  stiffness  A,  reduces 
7,/.  Reducing  the  inertia  /  reduces  y„. 

Investigating  the  ef  fect  of  shaf  t  speed  to  on  7  is  somewhat 
difficult  because  of  the  way  by  which  K,  and  D,  depend  on 
to  [see  Ref.  (#)].  However,  the  general  trend  can  be  obtained 
by  assuming  a  flexible  support  whose  stiffness  and  damping 
are  independent  of  the  shaft  speed,  e.g.  in  a  metal  bellows. 
For  this  case,  we  can  use  the  natural  frequency  of  the  system 
in  the  form 


to,,  = 

and  dehne  a  frequency  ratio 


.Adding  the  squares  of  the  lef  t-  and  right-hand  sides  of  [35], 
respectively,  gives 


1'^  + 

K  -I-  ~  n  'f 

4  ' 


[36] 


p  =  to/to„ 


and  a  damping  ratio 


«  = 


D* 

4/*io„ 


Multiplying  Eq.  [36]  by  [29]  gives  the  relative  transmissibility 


_ -f-  gf _ 

{K-if  -F  (o,  + 


(37) 


Ft|ualioris  (31.1],  [34],  and  [37]  give  the  relalicm  lietween  the 
relative  misalignment  7,  the  rotor  runout  7,.  the  initial  mis- 
.ilignment  of  the  stator  7,,.  and  the  dynamic  properties  of 
the  system. 


Replacing  the  dimensionless  parameters  in  [37]  with  their 
corresponding  dimensional  ones  for  the  case  D,  =  0,  A',  = 
A.^,  we  obtain 

_ _ 

(a;  -f  K*f,  -  /*co'')“  -F 

which,  by  using  the  aliove  definitions,  can  lie  arranged  in 
the  form 


DISCUSSION  OF  THE  RESULTS 


1: 

y> 


1(1 


[38] 


The  relative  misalignment  7  given  in  Fiq.  (34)  is  a  key 

factor  of  seal  performance.  This  misalignment  determines  Equation  [38]  is  well  dcK  iimented  in  many  text  IxMiks  on 

the  relative  pcrsilion  between  the  stator  and  rotor.  Parallel  vibration,  e.g.  Eq.  [3.2-4]  anti  Fig.  3.2-2  in  Rel.  (V),  pp  .')3- 
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'>  t  h  I  .o)  !m‘  dinru'ii\)at(Ml  will)  i<i  p  lo  dx' 

II i.i Mil) II  Ml  1 1  .lnsl))l^stl))li<  V 


(il,,,  ^ 


.11  llif  lre<|Ufni \  i alio 


P 


I _ 

(I  ‘-'t')' 


I  :<'->] 


|10| 


!•  I  oin  Ki|  I  tn|.  II  IS  » If. 1 1  ilial  i  lif  1 1  .iiisniissiluliu  jo  K<|.  |  .'<K  | 
li.is  .1  iii.iMiiiiiiii  oiiK  loi  •-  ill  wliiili  (asf  In  |;P.I) 

17  I  I'oi  t'  •  -.  ihf  II aiisniissiliilny  in  PW)  <I<h*s 

iioi  I  f.H  1 1  .111  f  XII  fin  II  III  .111(1  is  .ilwai  s  kss  than  11  nils  nanifl\ . 
7,,  7,  lor  an\  p  .\s  p  iiuifases  iiidflinilfk ,  (he  transniis- 

siliililv  III  IS.S]  appi oaflifs  nnil\  lor  anv  For  small  values 
ol  F.  ilif  ni.ixinuini  iiansmlssiliilits  in  |;I9|  lan  if a<  li  very 
lii^h  values  when  the  se.il  operales  near  its  lesonaiue 
()  I  !•  I  oni  the  .iho\ e  (lisi  iission.  it  tan  he  (oni  hided  that 
ihf  hfsi  (iioKf  ol  £  lot  l^efpin^  7,  niininiiiin  is  I/v2  <  i 
1.  whfif  £  -  I  (  01  responds  to  i  ritual  damping  whidi 
resulis  III  hisiesi  <lei.i\  ol  anv  disliirh.tiue. 

I  he  (  rilit.il  speed  to. ,  ol  la|.  [‘J  1 )  (an  he  expressed  in  let  nis 
ol  the  n.iiiii.il  ii('(|iifn<t  io„  gning 


(O 


I  III 


I  hiis.  the  ( I  ili(  .il  spft'd  IS  ,11  least  twi(  e  the  tt.titit  al  lre(|tieti(  v 
ol  the  ststetti  atid  (,iii  he  avoided  .iltogelher  hv  designing 
the  s(  .il  to  opei.ite  below  its  tiatiit  al  I  teipietti  \ .  This  (  an  he 
oht. lined  In  seleiiing  opliiiiinn  toning  that  iiiaximi/es  the 
llinil  him  silliness  [see  Rel  (7)|. 

As  w .IS  show  11  .iliov e  the  t  otor  t  iinoiit.  7,  h.is  .i  dii et  t  el  leet 
on  the  lel.itne  mis.ilignitient  7  wliu  h  (,ni  tanse  liiilitre  due 
to  lot.il  hue  (ont.ut,  Hetaiise  ol  maniil.K  tnritig  toleraiues. 

,1  (  ft  t. nil  I  inioni  is  mevitahle.  I  Inis.  ,111  nnpoilani  (|iieslioii 
III  in  in.  oiit.ii  ling  se.il  design  is  how  miu  h  rot 01  riinoni  t  ;iii 
tic  .illiiwed  hefoie  hue  (ont.ut  ottiits.  The  .niaivsis  pre- 
St  Iili  d  l.eie  W.is  h.ised  on  small  pel ttii hatioii  Ironi  the  po¬ 
sition  ol  p.n.illel  hues,  Kate  tonl.itl  represents  a  large  de- 
ti.iiion  lioiii  this  (ase  .ind.  henec.  may  he  iiiisnilahle  lor 
li.indling  hv  the  present  llieotv.  ^■el.  il  is  worthwhile  lo 
ex. inline  ihe  (oiuhiions  liii  hue  (oiinut  even  hv  the  tesnlls 
ol  the  sin. ill  pel  nil  h.iiion  .in.ilvsis  whi(  h  piovide  general 
lieiuls.  .md  .ilso  he(.nis('  when  .1  mole  el.ihoi.iie  an.ilvsis 
Ih'(  nines  .It  .nl.ihit  .  die  v  .iliilll  V  ol  the  I  elalivelv  simple  pres- 
i  lii  I  me  1  .III  he  et  .ihi. Ill'll. 

Kettiiinig  I  II.  III  Its  diiiiensionless  loiiii.  noting  that 
lot  sin. ill  |i(  I  till  h.iiion  .niv  .ixial  (list ni  h.nu  e  ol  the  stator 
dies  out  .diei  .1  while,  ihiis.  al  steady  stale  /.  -  0.  h  —  (■  ^ 
tve  h.ive 

//  I  +  7rt(osH  -t  (5  («  H.)  I  12| 

I'. lie  ( ont.ut  will  Inst  on  nr  al  H  ~  11  eilhei  al  the  inside  or 
onisule  i.idiiis  ol  the  seal  (le|K-n(ling  on  the  amonni  ol  Kin- 
nig  3  •  lie  ( (inditiotis  lor  hue  tonlad  are  then  Irom  |42] 
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I  7/h  -  (1 

lor  ((inlad  al  H  -  R„  and 

1  -  7  -I-  3(1  -«,)  =  0 


9 

I  lTi| 


[lUhl 


lor  ( onlat  t  M  R  ~  \ . 

Solving  K<js.  j4S]  lor  7  gives  two  dilleieni  values,  the 
smaller  ol  wliieli  delei  mines  the  appropriate  (iinlad  (on- 
(lilioii.  From  F.()s.  14!<|.  il  (an  he  seen  that  the  solution  lor 
7  IS  uletilual  in  the  vase  3  -  \/R,.  Foi  3  1'^/^..  lonlait 

(H  (  m  s  al  the  outside  1  adins  /{  =  1 .  w  hile  lor  3  >  1  ( oiitai  I 

odtirs  at  the  inside  radius  R  =  R,.  From  the  discussion 
lollowiiig  F',(|.  [20).  il  is  dear  that  for  the  case  I’..  >  P,  the 
opiimnm  dining  ^,^,1  =  'llR,(  I  -  R,)  is  larger  than  the  (  riliial 
one  3..  =  MR,.  Heiue.  lot  anv  reasonable  seal  design,  hue 
(ontad  may  (Hcnr  acdiiding  to  condition  |43a|  when  7  = 
\IR,.  From  Fap  (S4].  il  is  clear  that  7  reaches  its  niaxinuini 
v.ilue  at  T  ~  (1.  in  whidi  lase  7  =  7,,  +  7./.  Hence,  lor  3  '' 
3,,.  lace  (ontad  odiirs  when 


1"  '^'1  ^  ^ 

The  relative  tilt  7,,  is  related  to  the  rotor  runout  by  the 
relative  trunsmissibility  T  =  7,/7,  given  by  Eq.  [S7].  I'hus. 
F',(p  |44]  can  lx>  written  in  the  form 


7'(7.),-  +  7./ 


R, 


Hence,  the  critical  rotor  runout  causing  face  contact  is 


(7.),- 


I 


T 


[451 


where  7,/  is  given  by  F.q.  [51  a]  and  T  by  the  squ;ne  root  of 
K.q.  [57].  To  avoid  fat  e  (oniact.  the  critical  rotor  runout  (7,),, 
should  lie  as  large  as  |xissible.  From  [45],  it  is  clear  that  this 
will  Ire  achieved  by  revhtcing  7,/  and  T.  Hence,  the  same 
factors  that  reduce  the  relative  misalignment  7  increase  the 
critical  rotor  runout. 

In  diiM  lusion  ol  the  discussion,  a  numerical  example  will 
now  Ih-  presented  to  demonstrate  the  validity  of  the  small 
|H-rlurhaiion  assumption  in  practical  cases.  A  tvpital  seal 
with  the  lollovvitig  desigti  ptiratneters  atid  operating  (oti- 
ditioiiN  w.ts  seleded  .is  lollows: 


.Seal  outer  radiits,  »„  . 

Radius  ratio,  R,  . 

Design  dearaiu  e,  . 

Ftice  ta|>er  (cone  height )  . 

Stator  mass,  m*  . 

.Shaft  sjreed,  or  . 

Pressure  difleretitial,  p„  -  p,  . .  .  . 
Fluid  viscosity,  p.  (watet  at  liOT;) 

Supirort  stiffness,  K*  . 

Siip|Mirt  damping,  D*  . 


. 0.94  m 

. O.H 

.  10  - '•m 

. 2.5  X  It)  '>m 

.  "1  kR 

.  10  '  rad/s 

. 5  X  lo'Pa 

. 0.5  mPas 

. 5  X  lO’N/ni 

. .500  Ns/m 


1  he  Mippoil  stillness  and  damping  tocllitients  corre- 
s|><ind  to  a  Buna  N  secondary  seal  O-ring  at  the  given  shaft 
lie(|iiemv.  The  toriesponding  dimensionless  paiitmeteis 


lot  this  esaiiiple  .ire 

( ioning  angle.  (J . I2..'> 

Mass,  w  . 3.25  x  10  * 

Inertia,  I  . 1.63  x  lO"' 

Presstire  dilferential, /'„  -  /’,  . 0.26 

.Support  axial  stifTness,  A', :i  I  . 1.63  x  10  ' 

,Sup|)ort  angular  siiflness,  A . 6.14  x  10'^ 

Sitpport  axial  damping, /J.ei  . 9.76  x  10“  * 

Support  angulai  (lamping. /) . 4.88  x  10  '' 

Fluid  lilm  .ingulai  stillness.  A',  . 0.118  x  10  ' 

Fluid  lihii  aitgitlai  (lamping, /J/  . 0.83  x  10  ' 


Foi  this  ex.uiiple,  the  (alctilated  (iiti(al  speed  is(o,,  =  5.9(»>. 

1  he  liitio  7  ;  7  ,  =  0.06,  and  the  relative  transtnissihilitv  is 
y„y,  =  0.14.  Hem  e,  lot  a  rotor  rutiout  of  10  pm  at  the 
outet  dianielet  ol  the  seal,  which  is  a  lairlv  large  rtinoul 
( oi  responding  to  7,  =  1 ,  and  for  an  itiitial  tilt  of  the  stator 
ol  the  same  ordei  of  magnititde,  the  ntaxiinitm  relative 
misaligmnent  7  is  otilv  about  0.20.  It  should  he  noted  that, 
in  the  |)iesent  ex.itnple,  seventl  |iat.uneters  were  selected 
with  the  inletitioii  to  make  the  transtitissihility  oldrder  0.1. 
In  tn.itiv  pr.icti(.il  cases,  the  design  clettrance  and  shaft  six^ed 
will  he  sm.iller  atid  the  Ihtid  viscosity,  pressure  dif  ferential, 
support  .itid  Ihtid  lilm  stillness  atid  damping  will  be  larger. 
Iti  these  c.ises,  the  relative  transntissibility  wilLbe  even  less 
than  0.1  and  the  assumption  of' small  perturbation  more 
realistit  .  From  the  above  example,  it  can  be  seen  that  the 
(t'tti(,il  speed  (u,,  is  very  high.  Hence,  seal  stability  should 
tint  he  a  problem  iit  (otiventiottal  applications.  The  small 
pci  tut  hatiott  .tnahsis  predicts  the  stability  threshold  (|ui(e 
.IK  utiilelv  as  (ati  he  seen  Ironi  F.cjs.  123al  atid  (23b|.  Hence. 
F<i.  124|  lor  u).,  is  s.ilid  witliout  restrictioits  and  can  be  used 
to  examine  stable  operation  oi  any  noncontacting  seal. 

CONCLUDING  REMARKS 

1  he  dynamic  fiehavior  of  a  noncontacting  fttce  seal  was 
.m.il\/c(l  i.iking  into  account  manufaclui ing  and  assembly 
lolei  .im  es  in  the  form  of  rotor  ritnoiit  atid  stator  misalign- 
nicnl.  Complete  Ihtid  film  as  well  as  flexible  support  dy- 
nami(  cflcds  were  included.  The  etjuations  of  motion  in 
the  three  major  degrees  of  freedom  of  the  seal  were  pre¬ 
sented  and  solved  analytically.  The  solution  provides  a  sta- 
bilitv  threshold  in  terms  of  a  critical  speed  below  which 
stable  seal  operation  is  assured.  Steady-state  response  to 
lotoi  runout  in  the  presence  of  stator  initial  misalignment 
w.ts  lotind  and  expressed  in  terms  of  relative  transinissi- 
bilitv.  Flits  enables  ( ,il(  illation  of  the  lelalive  misalignment 
iH'iween  st.ilor  .ind  lolor  as  affected  by  the  manufaclurbig 
and  .issemblv  tolei.mces.  by  the  design  parameters,  and  by 
the  seal  operation  (ondilions.  It  was  shown  that,  in  the  pres- 
eiue  of  ImxIi  roloi  runout  and  stator  initial  misalignment, 
the  relative  misalignment  is  lime  dependent  and  varies  cvcli- 
cally  at  shaft  rotation  f  requency.  A  critical  rotor  runout  was 
found  related  to  various  design  parameters  and  operating 
((mditiotis.  This  ditical  value  can  give  the  designer  and 
manufacturer  ati  idea  of  the  magnitude  of  rotor  runout 


which  can  still  be  tolerated  without  causing  failure  bv  face 
contact. 

.'\  numeiit.il  example  w.is  piesented  lo  (l('nionsli .ite  the 
v.ilidilv  ol  the  sill. ill  pel  liii  b.ilioii  .issumplioii  in  pi.ulK.il 
.ipplualions.  Il  was  also  shown  ill. 11  slabilils  ihreshold  usu¬ 
ally  does  not  present  a  probletn  atid  its  predii  lion  hv  the 
small  |K-rlurbation  anah  sis  is  comparahle  with  ihisol  a  more 
elalxirate  numerical  one. 

Il  is  hoped  that  the  inloi  nialioii  proi  ided  here  will  dii  ed 
the  attention  of  seals  designers  and  users  to  the  various 
factors  affecting  nonconlai  ting  seals  perlorinance.  and  that, 
based  on  this  information,  (.ireful  experimentation  will  be 
carried  out  10  evaluate  the  iheoreiical  results. 
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APPENDIX  1 

Angular  Stiffness  and  Dannping  Coefficients  of  the 
Flexible  Support 

Usually  the  axial  stiffness  and  damping  coefficients  of  the 
flexible  support  are  measurable  and  known.  The  angular 
coelficienls  can  be  obtained  from  the  axial  ones.  Figure  8 
sliows  a  model  of  l lie  seal  in  wliicli  the  stator  has  a  displace- 
incnl  /  and  a  tilt  7,  alxiui  axis  v.  The  axial  displat  ement  ol 
a  |K>iiii  ti  ItHTHled  at  a  radius  r,  is 

/.II  ~  r,7,cosH  [46] 

where  r,  is  a  genet  al  raditil  location  of  the  flexible  siip|Kirl. 
and  6  is  given  by 

147) 


H  =  (fc  -  t|i 
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II 


1  he  axial  vehnilv  ol  point  H  is  theieloie 


y.ii  -  /.  +  /,  •yaosH  4-  i.y.tiisinO 


I'lH] 


By  definition,  we  find 


A.’  1 1 


d/Vfi 

ity 


Km, 


Du 

D\i  =  -  -  =  —  K 

ity,  2 


(5«1 


Kij 


m-. 

Hy. 


D^ 
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r'fOi 
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appendix  2 

The  Equations  of  Motion 

fii  Kef  .  (10),  a  kinematic  model  of  a  mechanical  seal  was 
pi  esented  and  the  dynamic  moments  for  the  case  of  a  flex¬ 
ibly  mounted  stator  with  antirotation  locks  were  obtained 
in  the  form: 


Due  to  the  motion  ol  |M)int  /<  an  axial  loi< e <lh'/  is  generaietl, 
its  value  is 

.//v  -  14‘1| 

where  (IKm\  and  .//fci  are  inlinitesimal  stiffness  and  damping 
loel  lii  ients.  res|)eclively.  Asstiming  uniform  tirettmferen- 
li.il  distribution  ol  the  tlvnamic  pro()erties  of  the  flexible 
support,  we  hate 


‘IK  1  ( 

K  ( (f/H 

L'tr 

I.WI 

'/Dll 

/>i:i'/H 

27T 

I  he  Miomenis  ai  ling  on  the  stator  due  to  die  flexible  snp- 
poi  I  ,ire  ihen: 

ahoiil  die  V  axis 


l\  =  /(7.  -  tit's.)  [50] 

7',  =  /(t|ii-y,  +  2tiry,) 

ft  was  alst)  shown  in  (10)  that  the  moment  is  of  order 
yf  and,  hence,  can  be  neglected. 

The  dynamic  moments  in  the  inertial  system  XFZ  are 
simply 

Tx  =  1%  [60] 

7'r  =  hy 

wliere  7,v  and  71  are  the  com|)onents  of  7,  along  the  inertial 
axes  ,Y  and  F.  respectiveb . 

The  angular  equations  of  motion  can  now  be  obtained  in 
the  rotating  system  xyz.  or  in  the  inertial  system  XYZ  by 
simply  equating  the  dynamic  moments  T  with  their  corre¬ 
sponding  applied  moments  Af.  Hence,  with  the  aid  of  Fig.  O, 
the  e<|ualions  of  motion  in  the  xyz  system  for  the  angular 
tiegrees  of  Ireedom  are 


.W  j  -  j  r.iosO'//'/ 

l-W] 

and  alMiiit  the  v  .ixis 

f’" 

Ml  -  J  r,sinH'/F/ 

(r.Hl 

Snhsiiiiiling  f  <p  1  l<i|.  .nid  1  liS|  dnoiigli  |.')  1 1  ,ind 
I'l'-’l  .111(1  I.M'I  gill's 

integraiing 

.i/i  -  -  ^lA'ii/ry,  +  /tiirry.) 

Af'j  ^  ^  />  i;iir’7,^ 

155) 

7(T  -  =  Af,x  +  Afpcosp  -  M/2sinp 

+  /Vf,Y,COSt|l  -I-  Af.,,COSp' 

/(>i|i7,  -Y  2tii7,)  =  Af„.  +  Afyisinp  +  Af/^cosp 
-  Af.visini|i  -  AfypSinp' 

In  the  inertial  .Yl'Z  system,  these  equations  are 
/•yv  -  M,\  f  Af/|(os<t>i  -  Af^',jsin4t|  +  Mx,  + 

1(12] 

lyy  =  M,y  +  Af/isin(t>i  -f  M/’icos^i  +  MYrSinoi/ 

where  M,x  and  M,r,  whicfi  are  the  Hexihle  support  induced 
iiurnienis  during  operation,  replace  the  moments  and 
M„.  These  new  introduced  moments  are  simply 
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A/,v  = 

-A.ii7v 

-  /Aii'Vv 

A/.i  - 

A,ii7i 

/f.ii7i 

rite  llui<l  liliii  iiidiiiL-d  inonirius  Af/j  and  Af/^  are  given  in 
K.qs,  (9a)  and  |9b|  in  terms  of  the  relative  misalignment  y. 
It  will  he  usclul  to  express  them  in  Kq.  [62]  in  terms  of  the 
angles  y\  and  yy.  From  Fig.  4,  we  can  see  that 


7\  =  "ycos*})!  +  7^osii)/ 
yy  =  7sin(t)i  +  •y,sin<i)/ 


[64] 


Dit lerentiatmg  with  respect  to  time,  we  have 

■yc<>s<t>i  -  'y<j>isin<t)|  =  y\  +  'y,sin<a/ 
■Vsin<J)|  +  -y<i>itos(t>i  -  -yi  -  y,cnso>l 


(651 


ly\  +  l)\]yx  +  A'ii7v  + 

1 66 1 

1 1<  i>sa>t  —  —  I  ixinti>t^  +  A, 117,1 
lyy  +  Diiyy  +  A'ii7>  -  -D/uyx 

[67] 

=  7r^A'iisin(i»Z  +  ^Z)/i|Cosa)Z 

F,quations  [66]  and  [67]  are  coupled  but  are  linear.  It  can 
also  be  seen  that  the  rotor  runout  7r  and  the  initial  mis¬ 
alignment  of  the  stator  7,,  lorm  forcing  functions  of  the 
dynamic  system. 

I'he  equation  of  motion  in  the  axial  degree  of  freedom 
is  simply 


Hence,  by  (oa),  (5b|,  [9a),  [9b],  [62],  and  [65]  F.qs.  [62]  can 

rearrangetl  in  the  form  tnZ  =  F/  [68] 
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t  he  (nmpli'le  iioiiliiteiii  equtiiiims  of  miilinn  tif  ii  fh’xihh  iiKiiinlrti 
sldliir  in  II  niiniiiiitni liiijf  loned-fme  iiiirhiiiiiiiil  si'iil  me  siilved 
iiiiineiiiiillx.  riie  '•nliilioii  itlilizei:  n  Iniiisieiil  il\iimiii<  miii!y\i.\  iiiul 
hikes  iiihi  iiiriiiiiil  nihn  (ixiril  ntiiiiiil  mill  iisseiiilils  liilermices  in 
the  liii  iii  iij  iiiiliiil  sliiloi  iiiisali^iiiiieiit.  (iiivihilinii  nf  the  fliiiil  film 
n  iilsii  m  run  nil'll  Ini.  .1  fimmlielrir  inveslif’itliiiii  is  fierf  iirmeil  mill 
the  efleil  nf  I'liimiis  ilrsi^ii  jimmiieleis  mill  iifHiiiliiiii  rnnililiiiiis 
nil  the  seal  ihniiiiiiis  is  jiieseiileil  mn!  ilisi  ii.sseil.  A  iriliciil  slutfl 
sfieeil  IS  liiiinil  nlnnv  u'hirh  the  seal  heniiiies  ilynmnirally  umtahle. 
.1  rnliral  rnlnr  niiniiil  is  fiiuinl  U'hirh.  if  exreeilril,  will  rame  seal 
Imliire  line  In  Inriil  fare  rubhinp;  riinlart.  A  aimparisiin  is  miule 
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between  the  nnmeriritl  irsiill.s  and  those  oj  a  simpler  anal'yliral 
sidnlinn.  It  is  fiiiind  that  the  mialytiral  sululnm  is  valid  fur  most 
prarliral  appliraliiins  of  merhaniral  seals. 

INTRODUCTION 

In  a  previous  pa|K’r  (/),  a  dynamic  analysis  based  on  small 
IK'i'lurbation  was  perrormed  lora  nonconlaclinirconed-race 
meclianical  seal.  The  analysis  i(M>k  into  account  rotor  axial 
runout  and  assembly  tolerances  in  the  form  of  initial  stator 
misalignment.  Both  stability  threshold  and  steady-state  re¬ 
sponse  of  the  seal  system  were  investigated.  A  critical  speed 
was  found  above  which  the  seal  becomes  unstable.  In  the 
stable  mode  of  o|>eration.  the  flexibly  mounted  stator  tracks 
the  rotor  with  a  certain  relative  misalignment  which  de- 
|>ends.  among  other  factors,  on  the  amount  of  rotor  runout 
and  assembly  tolerances.  For  large  runouts,  this  relative 
misalignment  can  lieconie  dangerously  high  and  cause  seal 
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(  =  (enn’i  line  <  learanrc 

t..  -  (Irsij^ii  (  U'.iraiue 

1)'  -  (l,iinpin^  MH'lfiueiil 
it  -  (liiiK’iisionIcss  (lamping 
A/  =  (linicnsioiilcss  axial  Uirre 

II  =  (liiiiensionk'ss  him  thickness,  liKt 
II,  =  lace  taper 

h  =  IcKal  him  thickness 
/  =  (litnensiotiless  moment  of  inertia 

K*  =  stillness  (oef  hcienl 
K  =  (litneiisionless  stiffness 

,\l  -  (limetisioiiless  moment 

III  =  (limeiisionli'ss  mass 

/’  '  (Imiciisioiilcss  piessiire,  pt.S 

p  =  pressure 

H  -  (limetisioiiless  ladial  ciMirdinate.  r/r„ 
>  ~  r.ulial  ( iHiiiliiiate 

■S  =  seal  paranielei .  (ijuo(r„/f,'„)‘  (I  -  II,)' 
I  =  lime 

/*  -  axial  (lisplacement 
/  =  dimensionless  displacement 


P*  =  lace  coning 
P  =  dimensionless  coning.  P*r,./f.'„ 

7*  =  relative  misalignnieiil 
7  =  dimensionless  niisalignmcni.  7*r„/fl„ 

7„  =  dimensionless  relative  iiiisaligiinieni  caused  hr  rotor  runout 

alone 

7,  =  dimensionless  rotor  runout 

7,  =  dimensionless  till  of  stator  (nutation) 

7„  =  stator  dimensionless  initial  misalignment 

7,;  =  stator  response  due  to  7,,  alone 

7.,  =  stator  response  due  to  7,  alone 

6  =  angular  coordinate 

p  =  visdisily 

ill  =  pretession 

(u  =  slialt  angular  veliKily 

Siilisi  1  i|Ms 

rr  -  ( ritnal 
f  -  llnid  film 
I  =  inner  radius 
m  =  mkl  radius 
»  =  outer  radius 

s  stator,  or  flexible  siip|Kirl 
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failure  due  to  excessive  leakage  or  local  lace  contact.  A 
critical  runout  was  found  in  (/)  which,  if  exceeded,  will 
result  in  lixal  lace  contact. 

I  lie  assumption  of  small  pertiirhation  f  rom  tfie  [Hisilion 
of  parallel  faces,  on  which  the  analysis  in  (/)  is  based,  leads 
to  a  litieari/ation  of  the  ec|uations  of  motion  of  the  seal  ring. 
I  his  enaliles  c  losed-loi  tn  atialytical  soltition  and  yields  very 
useful  analytical  expressions  for  seal  civnamic  fjerformance. 
Howevei.  the  small  |)et turhation  assumption  may  lie  re¬ 
st  ric  tive  in  c  ases  wheie  the  laces  are  fat  f  rom  lieing  parallel. 
I  lie  case  of  c  ritic .il  rotoi  l  unoiit.  lot  example,  may  lie  in- 
,ippiopi  i.iic  lor  li.mclliiig  l>\  tlie  small  pei  lurhation  analysis 
smie  it  coiiesponds  to  high  lel.itiye  niis,ilignmeni  lielween 
lhese,il  l.ices  .\iioihei  lestric  lion  of  tfie  small  |ierlurl)alion 
assuiiiptioti  is  the  tiegfecl  of  c  ayitation  and  nonlinear  el  feds 
on  the  dvtiamic  pioperlies  of  the  fltiiil  film.  Meiice,  high¬ 
speed  locs-ptessurc-  cases  may  he  ill  lio.tled  fiv  the  simple 
analysis.  In  otcler  to  cover  llie  full  range  of  seal  operation 
conditions,  .tticl  to  olii.tiii  a  complete  pic  lure  of  its  dynamic 
lietiayior.  the  small  jjerturhaticin  restriction  has  to  be  re¬ 
laxed.  and  all  the  nonlinear  ef  fects  liave  to  lie  accounted 
led  .  1  liis  i  ei|uii  es  .1  miiiiei  ic  .il  ir.iiisieiit  solution  of  itie  com¬ 
plete  d\ti. unit  et|it.tiiotis  of  the  seal.  It  ts  the  pitrpose  of  iliis 
lepotl  to  pic'setil  the  tesitlls  of  site  h  ati  elfotl.  Itt  addilitm. 
a  c  ottipat  isoti  will  Ik'  itiacle  with  tfie  testtlts  olilaitied  by  the 
small  iM't'titt liatiott  .ttialysis  of  Kef.  (/).  I'liis  will  determine 
the  range  of  valichty  of  tfie  simpler  and  more  conyenienl 
to  use  anahsis  of  (/). 


a  few  milliradians,  and,  hence,  can  fie  treated  as  veciors. 
•Since  y*  is  fixed  in  space  and  yf  is  rotating  at  the  shall  speed 
CO,  the  resultant  vector  yf  will  possess  a  time-yarying  speed, 
ifi.  In  Ref,  (/),  il  was  found  that  the  vector  y*  can  be  ex¬ 
pressed  in  the  dimensionless  form  (see  Nomenclature  for 
y*  and  -y): 

+  y.,  ni 

where  7,/  is  hxed  and  is  the  response  to  alone,  and  7,, 
is  rotating  and  is  the  res(X)nse  to  y,  alone.  As  wifi  lie  shown 
later,  a  similar  presentation  of  7.  is  possible  in  tfie  nonlinear 
case,  tcM).  The  relative  misalignment,  7,  between  the  stator 
and  rotor  is  also  a  rotating  vector  given  by  the  vector  sub¬ 
traction 

7  =  7.  -  [2] 

Figure  2  shows  the  relative  position  between  the  seal  com¬ 
ponents  and  Fig.  S  gives  a  vector  representation  of  the  var¬ 
ious  tilts  of  the  seal  system.  The  vector  y„  in  Fig.  3  is  the 
relative  misalignment  7  in  the  special  case  7,/  =  0  and  is 
simply; 

7..  =  7.-  -  7r  [3] 

Fhe  eejuations  of  motion  of  the  Hexibly  mounted  stator 
are  (see  Ref.  (/)]: 


BACKGROUND  AND  METHOD  OF  SOLUTION 

rhe  seal  model  is  shown  in  Fig.  I.  Il  consists  of  a  seal 
se.il  (rotor)  that  is  rigidly  mounted  to  ilie  rotating  rigid 
sli.dl.  .mil  a  flexilib  supported  seal  ring  (stator).  I'he  rotor 
h.is  an  axial  runout  that  can  lie  represented  by  a  lilt  yf  of 
Us  plane  with  respec  1  to  the  axis  of  shall  rotation,  bimilariv. 
I  lie  stator  mav  have,  prior  to  hnal  attachment  to  the  rotor 
(see  (/)],  an  initial  misalignment  7*  with  respect  to  the  axis 
of  shaft  rotation.  ,-\t  rest,  and  with  zero  pre.ssure  differ- 
ciilial.  the  stator  is  pressed  against  the  rotor  by  the  sup¬ 
porting  springs  alone.  This  forces  the  stator  into  the  same 
nil.  y*.  as  that  of  the  roloi  .  During  ci|>eralion,  the  mating 
l.ii  es  separ.iie  and  the  stator  assumes  a  new  till  7.*.  l  itis  till 
IS  a  result  of  the  combined  effects  of  Ixilh  7,*  and  7*.  Fhe 
till  angles  7*.  7*,  and  7,'  are  all  very  small,  of  the  order  of 


/(7,  -  «i<S.)  =  Afx 

[4] 

/(il/y,  +  2ii)7,)  =  My 

[5] 

mZ  =  F/ 

lb] 

where  M,  and  M,  are  the  moments  acting  on  the  stator  afxiut 
axes  X  and  y,  respectively,  of  a  reference  system  xyi  which 
rotates  at  a  rate  ij;  in  an  inertial  system  XVZ  [see  Figs.  2(a) 
and  3).  I  he  axis  x  of  the  rotating  system  coincides  with  the 
vector  7,  and  is  at  an  angle  i|>  with  respect  ni  the  inertial 
axis  X. 

The  moments  an<l  Af,  as  well  as  the  axial  force  /v 
consist  of  (ontribiiiions  f  rom  Ixith  the  flexible  support  and 
the  fluid  lilm.  I  he  supjxirl  contribution  to  the  moments 
and  ftirces  is  (/): 

=  A.',  [7„c<>si|i  -I-  7,<  t>s(i|>  -  gg)  -  7,)  -  D,y,  [7] 
Af.,  =  -A’,[7„simji  +  7,sin(«li  -  <iX)]  -  D.tjry,  (8) 

F.Z  =  -k:z  -  d:z  [rj 

where  A','  and  D[  are  axial  stiffness  and  damping  coefficients, 
respectively,  of  the  support .  T  hese  coefficients  include  the 
dynamic  properties  of  the  supporting  springs  and  the  sec¬ 
ondary  seal  which  can  be  a  metal  or  elastomer  bellows  or 
elastomer  ring.  The  dynamic  pniperties  of  an  elastomer 


NoiiliiK'ar  Dvtiainu  Analysis  ot  NonciuUarliii)'  (aincd-Kad'  Met  lianical  Seals 


h 


Fig.  2(a>— SmI  modal  and  coordinala  ayatama 


Fig.  3— Vactor  rapiaaantatlon  of  tIHa  of  tha  aaal  eomponanta 


se<<in<lary  seal  t  an  lie  obtained  experimentally,  e.g.  (2).  From 
the  known  eoellicienls  and  I>1,  one  can  calculate  the 
angular  slil loess  and  damping  coeilicients  K,  and  1),  (/). 
The  terms  y,  and  y,,  in  Kqs.  [7]  and  |8]  are  due  to  the 
assembly  tolerances  (y,,  is  arbitrarily  assumed  to  coincide 
with  the  negative  X  direction). 

The  Huid  Him  contribution  to  Af„  and  F/  is  obtained 
by  Hrst  calculating  the  pressure  distribution  in  the  sealing 
gap  and  then  integrating  numerically  over  the  face  area. 
Expressions  for  the  various  fluid  film  pressure  components 
were  given  elsewhere,  e.g.  (5),  and  are  not  repeated  here. 
It  should  be  noted,  however,  that  cavitation  is  accounted 
for  in  the  present  analysis  by  deleting  any  negative  pressure 
Itefore  integration.  In  addition,  curvature  effect  due  to  the 
change  in  r  is  alst)  accounted  for  in  the  numerical  integra¬ 
tion. 

I'he  dimensionless  local  film  thickness  is  given  by 

H  =  I  +  E«cose  +  6(«  -  ft,)  [10] 


where 


and 


e 


y 

l+Z 


a 


P 

i+z 


•Since  y  and  Z  are  time  dependent,  H  is  time  dependent, 

t<H>. 

The  relationship  between  the  various  tilt  angles  is  (see 
Fig.  ,S) 


SECTION:  A-A 


Fig.  a(b)— NsMfv*  poaMen  hafwaan  mM  ring  and  aaal  aaal 


■y  =  [7?  +  7r  -  2y.y,cos(t|i  -  <•)<)]''*  [II] 

Hence,  starting  with  some  arbitrary  values  for  the  initial 
conditions  y„  i|i,  Z  anti  their  corresponding  time  derivatives, 
the  him  thickness  H  can  l>e  calculated  by  Eqs.  ( lOj  and  [  1 1 1. 
Next,  ilie  pressure  lielil  anti  the  resulting  mttmenis  ,Vf, 
and  the  forte  /•>  are  i  alt  iilated.  I'he  time  seiftiid  derivatives 


D-4 
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(>(  llu-  sl.itoi  ol  li  cediitii:  y,.  ij).  and  /.  .iir  ihfii  tiiiind 

Ironi  K(|s.  I  t|  1(1  |(i|  III  lilt-  tonii: 


ek-iiifiil.  'Mif  fX|)rfssioii  lor  K*  and  D*  arc  (6): 


y  +  <l»  y, 

[12] 

Af,  .  \ 

T  - 

[131 

,  /•/ 

/  =  — 

[14] 

K*  =  Mpo  -  p.)  p-  m.  -  1) 

Df  =  12iT^ip(l  -  R.fRJO. 


1  -  R,  ’ 

4  +  23(1  -  R,) 


These  values  are  inlegraterl  with  respect  to  lime  and  new 
values  are  Tound  hit  y„  y„  ij),  X  and  Z.  The  niinimiim 
him  thickness  is  <  he<  ked  at  each  time  step  and  the  proce¬ 
dure  repeats  itsell  until  liKal  face  contact  occurs,  or  until 
the  minimum  him  thickness  shows  steadv-siate  liehavior  over 
a  lew  suitessive  sli.ill  revolutions.  The  (omptiter  program 
used  is  hasicallv  sittiilar  to  that  descrilH’d  itt  (•/)  with  a  Tew 
niodihi.iiiinis.  The  lime  inlegration  inenlioned  alxive  was 
pet  lot  I  lied  h\  the  M 1  l.N  T.  inlegration  nielluHl  ol  the  C.SMl’ 
(^)  on  an  IhM  ;17(I/I(ih  tomputei.  The  average  ttmiptiler 
lime  reiiuired  was  ahoul  1  s  pet  each  revolution  oTthe  rotor, 
and  it  took  at  least  six  revolutions  to  determine  the  stability 
status  of  the  seal  svstem  Tor  a  given  set  of  operation  con¬ 
ditions  and  design  parameters. 

RESULTS  AND  DISCUSSION 

The  computer  program  enables  one  to  de.scrilie  the  com¬ 
plete  dvnatiiic  behas  ioi  ol  the  seal.  This  includes  both  stator 
tracking  oT  rotor  runout  in  the  stable  mode  of  operation 
and  dvnamic  slabiliu  threshold  oT  the  seal.  As  was  shown 
in  (/),  the  axial  runout  y,  and  the  initial  misalignment  oT 
die  stator.  y,„  Torm  torcing  Tunclions  and,  hence,  af  fect  only 
tile  dvnamic  resjxmse  at  stable  operation  but  not  the  dy- 
ii.imic  stability.  Yet,  it  is  interesting  to  compare  the  stability 
threshold  obtained  Irom  the  numeriial  .inalysis  with  that 
iihlained  from  the  a|)proximale  analytical  solution  of  Ref.  (/). 

.\  i  rilical  speetl  to,,  was  Touiul  in  Ref.  (/)  almve  which  the 
se.il  Iteiomes  unstable  (stability  threshold).  The  analytical 
expression  foi  w,,  in  (/)  is 

=^(‘  ^^)(' 

In  Eq.  ( I.*)]  A'*  and  W,*  are  the  angular  stiffness  and  damp¬ 
ing  coefficients,  respectively,  of  the  support  related  to  the 
axial  coefficients,  K*'  and  Df  by  [see  (/)]; 


K*  =  a:- 


<)i[l  -r3(l  -«,)]  -  2  - 


3(1 


D*  =  nr 


where  r,  is  the  radial  ItKation  of  the  Ilexibic  sup|ioriing 


_  -.O  -2  4-  3(1 -A,) 

~  3 

The  preiession.  ol  the  st.itor  at  the  stability  threshold 
was  found  in  (f)  in  the  torni 


A  stability  threshold  search  was  performed,  using  the  pres¬ 
ent  computer  program  for  60  different  cases  to  evaluate 
the  accuracy  of  the  small  perturbation  analysis  in  predicting 
seal  stability.  The  range  of  radius  ratios  covered  was  be¬ 
tween  R,  -  0.7  to  R,  =  0.98.  At  each  radius  ratio,  several 
combinations  of  K*IK*  and  D*/D*  values  were  tested.  These 
combinations  consisted  of  the  values  0.2,  0.4,  0.6,  and  0.8 
for  K*/K*,  and  the  values  0,  0,2,  and  0.4  for  D*ID*.  A 
relatively  low  pressure  differential,  having  the  dimension¬ 
less  value  P„  -P,  =  0.25,  was  selected  to  facilitate  unstable 
o|>eration.  The  face  coning  3  at  each  radius  ratio  was  the 
optimal  one  for  maximum  K*.  and  is  given  by  [see  (6)]: 


R,(l-R.) 

The  numerical  analysis  gave  critical  speed  values  that  were 
always  slightly  less  than  those  given  by  Eq.  [15].  The  dif¬ 
ferences  ranged  from  4  percent  at  R,  =  0.7  to  only  1  per¬ 
cent  at  R,  =  0.98.  The  results  for  tji.r  were  identical  to  those 
of  Eq.  [16].  Hence,  it  can  be  concluded  that  the  small  per¬ 
turbation  analysis  of  Ref.  (/)  predicts  quite  accurately  the 
stability  threshold  of  the  seal  and  Eq.  [  1 5]  can  be  used  safely 
to  calculate  the  critical  speed. 

The  effect  of  rotor  axial  runout  and  stator  initial  mis¬ 
alignment  on  the  dynamic  behavior  is  fully  expressed  in  the 
stable  mode  of  seal  operation.  The  computer  program  was 
used  to  calculate  the  minimum  film  thickness  (A/C„),„,,  the 
relative  misalignment  (and,  hence  the  transmissibility 
yjy,),  and  (he  axial  displacement  Z.  of  the  flexibly  mounted 
stator.  'This  was  done  for  various  values  of  nrtor  axial  run¬ 
out,  7„  and  ft>r  =  0.  A  very  important  result  of  these 


It  • '  ^  ^  . 
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I  alt  iil.ilioiis  Is  (lie  irilual  rtiloi  niiitiul  {y,),,  whitli  i.uises 
liH.il  lat e  I onlat  I  ii.iiiiclv.  (/</^.„)„„„  0  Knowlt-tl^i' iii  llie 

tiiliial  M>ii)i  runout  is  very  Ir'I|>IuI  in  setliii^  reasonable 
inannlaetnrin^  loleraiues  and  in  avoiding  seal  lailnre.  Ki- 
nallv.  die  ellett  oi  assembly  misalignment,  y,„  on  ilie  dy- 
naniu  Irebavior  was  also  examined.  I  bis  misalignment  pre¬ 
vents  tracking  with  unilorm  nutation  and  precession.  Hence, 
the  relative  misalignment  lietween  stator  and  rotor  as  well 
as  the  axial  clearance  liecome  time  depetident  and  vary 
( vchcallv  at  the  shaft  f  re<|nenrv. 

.'\ltogether  theie  are  ten  dillerent  parameters  whith  af- 
let  t  the  seal  tbiiamits.  These  .ire:  the  seal  outer  and  inner 
r.idii.  stator  mass.  Ilnid  \  isiositv,  piessnre  dil  leremial  .icross 
the  seal,  designetl  axial  tle.n.nue.  I.ite  tonnig.  sh.d/speerl. 
.md  siiftness  and  damping  of  the  flexible  snppoil.  .\  para- 
metrit  inxestig.ition  whith  lakes  imo  attonm  variations  in 
.ill  these  parameters  wtinid  re<|nire  .in  enoi  inons  computer 
time.  Hence,  to  overcome  itiis  problem  wilhtiut  losing  ttK> 
imu  h  inlormation,  a  ‘single  pertnrfiaiion"  method  was  uii- 
li/ed.  Bv  this  method,  a  reference  t  ase  is  set  np  for  the  ten 
parameters.  Kacti  of  the  more  import.mt  parameters  is  per¬ 
turbed  one  at  a  lime  and  then  regains  its  refereiue  value. 
Solutions  ate  obtained  lot  each  single  peiturbalion  result¬ 
ing  in  reasonable  amount  of  inform.ilion  from  a  relalivelv 
small  numiter  ol  v.iriations.  Mans  cases  wfie  examined  by 
the  analysis  of  (/)  jnior  to  selecting  the  referetue  tase.  In 
most  practital  tases,  the  iransmissibililv  was  fbutid  to  Ire 
yjy,  I.  Hence,  most  prat  lit  al  cases  t  tmiply  with  the  def¬ 
inition  of  small  ))eruirbation  lor  which  the  results  tif  Ref.  (/) 
shtiuld  Ih’  valitl  This  was.  indeetl.  verified  bv  solving  some 
of  these  t.ises  nnmeiitallv  .iiitl  tomp.u'ing  the  results  with 
those  of  Ref.  i/l.  I  he  loiielalion  w.is  always  within  a  lew 
pen  em. 

In  Older  to  deieimint'  the  i.nige  ol  v.ihdilv  of  the  small 
pel  inrb.iiion  .ui.ilvsis  ol  ( / ),  t  .ises  w  hit  h  were  t  learly  tun  til 
die  I.nige  ol  small  perinrbalion  were  singled  out  anti  se- 
leitetl  loi  the  reference  se.il.  These  tases  are  t  harat  leri/.etl 
b\  rel.itiveb  large  mass,  large  axial  t  learant  e,  very  high  shaf  t 
s|K-etl.  low  pressure  differenlial.  .mtl  low  fluid  vist  tisitv.  The 
se.il  seletletl  tor  the  reference  case  tit  the  parametric  in- 
\esiig.iiion  is  a  nietliiim  si/e  seal  with  BO  mm  outer  diameter, 
.ind  .1  O  H  i.itlins  ratio,  sealing  water  .it  lem|H’raluie  of  ()0°('. 
.mtl  piessnre  ol  O  .'f  MI’a.  I'he  staloi  mass  is  I  kg  and  the 


tlesignetl  sealing  g.ip  is  H  )ini.  The  stillness  anti  ilampmg 
of  die  snp|Mirl  aie  ivpit.il  lor  Buna  N  O-ring  (2).  An  o|>- 
tiimiin  ttuiing  with  20-p.ni  taper  was  selet  ted.  Table  1  pre¬ 
sents  the  referetue  t  ase  and  the  jierturbations  in  the  more 
important  paranieteis.  I  he  stiffness  anti  damping  of  the 
Ilexilile  su|>|M>rt  were  int  reasetl  oi  det  reasetl  simidlane- 
tiuslv.  I  he  high  values  t  tirres|>tind  tti  X’ittin  75  O-ring  while 
the  low  values  lepreseni  a  metal-lielltiws  sectindary  seal. 

The  results  lor  the  case  y„  =  0  are  shown  in  Tigs.  -I  thrtmgh 
K.  Tarts  (a),  (b),  and  (t )  til  each  figure  present  the  mininmni 
film  thickness  the  relative  misalignment  y„.  and 

the  axial  tlisplacement  /.  res|x-ttively.  as  lunttitins  of  the 
rottir  runout  y,.  \  t  iiniparison  is  made  lietween  the  present 
nniiierical  stilmitin  anti  the  analytical  stilution  of  Kef.  (/). 
Trtini  the  part  (a)  of  eat  h  ol  the  Tigs.  4  through  B.  the  t  ritital 
rtittir  riintiul  t  an  Im'  obtained  frtini  the  intersectitin  of 
the  t  tirves  with  the  line  f/i/f.',,),,,,,  =  0.  It  is  quite  dif  ficult  to 
tibtain  exactly  this  intersectitin  frtini  the  numerical  stilmitin 
tiue  tti  the  very  large  ctimputer  lime  required  in  cases  tif 
large  runouts.  However,  the  intersection  can  be  obtained 
from  the  available  numerical  results  by  exira|Milalitin.  It 
shtiuld  also  lie  noletl  that  as  the  rotor  runtiut  approaches 
its  critical  value,  the  axial  tlisplacement  7.  Iiectimes  very  large 
[see  part  (t )  of  the  figures).  Hence,  in  many  cases,  the  seal 
may  fail  lu’cause  of  an  ext  essive  leakage  much  lieliire  Itual 
face  ctintact  tRcurs.  It  is  also  more  practical  tti  set  a  limit 
for  the  |>ermissihle  minimum  Him  thickness  at,  say,  f/i/ 
=  (1.25,  and  dehne  the  critical  rotor  runout  as  the  value 
of  "y,  ctirresptinding  to  that  limit.  This  will  prevent  the  po.s- 
sibiliiy  of  Itical  facecomaci  and  will  provide  a  facttir  of  safety 
and  a  intire  practical  value  ftir  the  critical  rotor  runout. 
There  were  no  tliflit  iiliies  in  obtaining  numerical  solutitins 
lor  t  ases  tif  (li/C,,),,,,,,  »  <»  '25. 

Figure  4  presents  die  elfet  I  of  the  design  t  learant:e  f.„. 
Krtim  Tig.  4(a).  it  can  be  seen  that  the  critical  rotor  runtiut 
obtained  from  the  numerical  solution  is  always  larger  than 
that  predicted  analytically.  This  is  true  fur  all  the  other  cases 
as  well  (see  part  (a)  tif  Figs.  5-8].  Hence,  the  analytical 
prediction  may  serve  as  an  upper  limit  ftir  the  critical  runout 
which,  if  not  exceeded,  assures  noncontacting  tiperatititi. 
The  analytical  (y,),,  values  predicted  for  design  clearances 
(:„  tif  12  pm,  8  pm,  and  4  pm  are:  1.745,  3.310.  and  8.200. 
ies|>ectively  (Fig.  4(a)l.  The  axial  runout  at  the  outer  radius 


*  m"*  / 


Cl 


. -f.  4 

^  ' 


r*n 


Noiilinc.ii  I)Mi.iiiii>  Aii.ilvsis  <>l  N'liiuoiit.u  linK  (liiiifd-I  .Kt-  Mt-i  li.iiiii  al  Seals 


=  0.25,  (he  analytical  results  Ibr  the  corresponding 
y,  are  from  20  to  30  percent  smaller  than  the  numerical 
lesults  depending  on  the  value  of  p„  -  fi,.  The  lietter  ac- 
(  iirat  V  Kit  responds  to  (he  highei  pressure.  From  Fig.  (i(h), 
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Fig.  5<a)— Minimum  film  thicknMS,  (fi/Co)mn.  v<  rotor  runout,  y„  and  aliaft 
spaed,  ui. 
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it  can  be  seen  that  increasing  the  pressure  dilierential  re¬ 
duces  the  iransmissibility,  and  hence,  reduces  the  danger 
oi'liK'al  lace  contact.  F'igure  ()(c)  shows  that  higher  pressure 
dilferential  redutes  the  axial  displacement  /.  Fliis  is  proh- 
al>ly  due  to  the  suppression  ol  cavitation  eliect  which  in¬ 
duces  liydrtKlynamic  axial  lorce  tending  to  open  the  sealing 
K^p. 

F'igure  7  presents  the  effect  of  various  flexible  supports. 
Case  I  corresponds  to  a  metal  l)ellows,  case  II  to  a  Buna  N 
O-ring  secondary  seal,  and  case  Ill  to  a  Viton  73  f)-ring 
secondary  seal.  From  Fig.  7(a),  it  can  be  seen  that  the  Buna 
N  which  provides  the  medium  stiffness  and  damping  is  the 
best  choice  at  the  ojrerating  speed  of  the  reference  case. 
The  Viton  73  gives  the  best  correlation  between  the  ana¬ 
lytical  and  numerical  results  (about  18  percent  at 
(/i/C„)„,„  =  0.25).  From  Fig.  7(b),  we  see  that  the  iransmis- 
sihility  y.Jy-  *('  almost  unaffected  by  the  variations  in  the 
flexible  support.  A  similar  conclusion  can  be  drawn  from 
Fig.  7(c)  regarding  the  effect  on  the  axial  displacement. 

Figure  8  presents  the  effect  of  lace  taper,  H,.  The  taper 
is  related  to  the  lace  coning  by 

//,  =  3*(r„  ~  r,)  [18] 

Hence,  by  F.qs.  [18]  and  [17],  the  optimum  face  taper  is 


(-]  -- 


and  in  our  reference  case  for  R,  =  0.8  and  C„  =  8  p.m  (see 
Table  I)  the  optimum  ta|)cr  is  H,  =  20p.m. 

From  Fig.  8(a),  we  see  that  as  the  taper  decreases  the 
ininiiniim  film  thickness  inc  teases;  also,  the  correlation  be¬ 
tween  the  analytical  and  numerical  solutions  improves.  Fig¬ 
ure  8(b)  shows  that  reducing  the  taper  reduces  the  trans- 
inissibility.  and  Fig.  8(c)  shows  that  the  smallest  axial 
displacement  correspontls  to  the  smallest  taper.  These  re¬ 
sults  can  be  attributed  to  the  increase  in  the  fluid  film  damp¬ 
ing,  1)/,  as  the  face  coning  decreases  (6).  From  Fig.  8,  it  seems 
as  il  the  smallest  taper  of  10  pm  is  a  better  choice  than  the 
opliimim  coiling  whicii  in  our  case  is  20pm.  This,  however, 
is  iiiisifacliiig  lx.'caiise  ol  the  effect  ol  the  coning  on  the 
stability  iliresliold  and  on  the  transition  |)eriixl  lietore  steads 
Slate  is  reached.  As  was  slated  before,  small  coning  results 
in  high  fluid  lilm  damping  and,  hence,  reduces  the  critical 
s|jeed  Cl),,  (see  Kq.  [13]).  For  the  reference  case  of  Table  1, 
the  ratio  of  critical  s|>eefl  to  operating  speed,  lo, ,/(!),  is  2.5 
for  H,  =  40  pm,  2.51  for  H,  =  20  pm,  and  2.29  for 
W,  =  10  pm.  Kediicing  llie  la|>er  to  H,  =  3  pm,  for  ex¬ 
ample,  results  in  ci>,,/ci>  =  1.20  and  the  .seal  approaches  its 
stability  threshold.  Further  reduction  of  the  la|>er  would 
cau.se  seal  failure  due  to  dynamic  instability. 

Another  disadvantage  of  small  taper,  which  in  practical 
seals,  is  associated  with  high  fluid  film  damping  exceeding 
the  rrilical  value,  is  the  slow  decay  of  any  disturbance  of 
the  dynamic  system.  In  order  to  obtain  the  steady-stale  re¬ 
sponse  shown  in  Fig.  8  at  -y,  =  1.5,  for  example,  16  shaft 
revohnions  were  rei|iiirc‘d  for  the  case  H,  =  10  pm  com- 
|Kire<l  to  only  8  and  5  revolutions  for  the  20-pm  and 
40-pni  ia|H'i's,  res|H‘ciively.  Increasing  the  rotor  riinoui  ni 
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Fig.  6<a) — Minimum  film  thicknau,  (h/CoJmm.  va  rotor  rutKMit,  y„  and 
praaaure  differantlal,  Po-pi. 
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Fig.  6(b>— Relative  tilt,  y,  va  rotor  runout,  y,.  and  praaaure  differential, 

Po-Pi. 
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Fig.  a<c>— Axial  dlaplacemant,  Z,  va  rotor  runout,  y„  and  preaaure  dlf- 
farantlal,  Pa-pi. 


y,  =  ■>.<),  (<»r  example,  rerpiired  1 1  shaft  revolutions  to  reath 
steady  state  at  H,  -  20  |im,  whereas  at  H„  -  10  |im,  even 
.50  revolutions  were  not  enough.  Hence,  it  can  he  concluded 
that  the  optimum  coning  would  he  a  lieiter  choice. 

In  general.  Figs.  4  to  8  show  that  reducing  the  rotor 
runout  7r  increases  the  minimum  film  thickness  and  reduces 
both  the  transmissibility  and  the  axial  displacement.  This 
results  in  a  seal  that  operates  with  its  faces  closer  to  parallel 


|Hisilion  and  reduces  die  tlanger  of  f  ailure  due  lo  loc  al  lace 
contact  or  excessive  leakage.  .\s  the  rotor  runoui  inc Teases, 
the  axial  displacemeni  increases,  too.  Siiue  die  leakage  de- 
|>ends  on  f.'*.  and  f.'  =  (,„(  1  +/),  an  increase  in  /  has  a  verv 
strong  effect  oti  ifie  leakage,  fl  seems  that  mani  seals  mai 
fail  tlite  to  high  leakage  much  befoi  e  f at  e  i  uhhing  as  a  result 
of  t  ritical  lolot  i  iinout  oci  tirs. 

/Xnothet  important  finding  from  Figs.  4  to  8  is  the  fait 
that  the  .uialvtiial  prediition  for  die  minimum  film  tliiik- 
iiess  ihJ(. IS  ahvavs  belou  the  numei  kal  result  anil.  Iieiiic, 
is  on  the  s.ile  siile.  I  bis  .ilso  means  that  the  i  ritital  rotor 
iiinoiit  found  in  (/)  in  the  loiiii 
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(ivliete  /'  IS  die  iraiisniissibihiv  yjy,t  is  a  sale  limit  lor  tlesign 
pur|M>ses 

III  order  to  examine  the  elfeil  of  assemblv  tolerances  on 
the  tiynamic  behavior  of  the  seal,  an  initial  stator  misalign¬ 
ment  of  7„  =  5  was  selectetl.  I'liis  misalignment  corre- 
s|>onds  lo  O.li  mrad  in  the  lelerenie  case  of  Fable  1.  With 
this  misalignment  added  lo  the  reference  case,  solutions 
were  ohiainetl  al  rotor  runoui  values,  y,.  of  ;  0.5,  0.75.  and 
1.5  (torresponding  to  0.00,  0.15,  and  0.5  mrad.  respec¬ 
tively).  File  anaivsis  ol  Rel.  (/)  was  u.sed  lo  calculate  the 
stator  resjMin.se  7,/  ami  it  was  found  as  7,/  =  0.120  corre- 
s|>onding  to  0.020  mrad  in  our  reference  seal.  The  critical 
rotor  runoui  (7,),,  in  the  presence  of  the  stator  initial  inis- 
aligninent  was  calculated  [using  (/)]  and  found  to  lie  10 
|)ercent  less  than  the  critical  value  at  7,,  =  0.  Hence,  initial 
stator  misaligiiment  redut  es  the  critical  rotor  runoui  as  can. 
indeed,  lie  seen  from  Ki|.  [20|. 

Figures  0  through  1 1  jireseiii  the  steady-slate  res|M)nsc 
of  the  flexibly  mounted  stator  in  the  presence  of  Ixith  rotor 
runout  7,  ami  initial  stator  misalignment,  7,,.  In  Fig.  9,  the 
axial  motion  Z  is  shown,  while  Fig.  10  presents  stator  nu¬ 
tation  7„  and  Fig.  1 1  shows  the  variation  in  the  phase  angle 
t|i  -  u>/  (see  Fig.  5)  between  the  tilt  vector  of  the  staloi  7..  and 
the  till  vector  of  the  rotor,  7,.  Fhe  time  base  in  these  figures 
is  expressetl  in  shaft  revolutions.  From  the  three  figures,  it 
can  Ik’  seen  that  the  stator  resptinse  consists  of  cyclic  vari¬ 
ations  in  /.  7„  and  the  phase  angle,  and  that  the  fre(|ucncy 
of  these  variations  is  etpial  to  the  frequency  of  the  shaft 
rotation.  .Such  Itehavior  was  already  predicted  in  (/)  and 
afso  olxserved  ex|)erimentally  (7),  (^).  It  should  be  noted 
here  that  while  the  small  |>erturbation  analysis  (7)  can  pre¬ 
dict  the  variations  in  7,  and  the  phase  angle  between  7,  and 
7,.  the  lime  variation  in  Z,  which  is  due  to  the  coupling 
between  the  axial  and  angular  degrees  of  freedom,  can  be 
predicted  only  by  the  numerical  solution  which  takes  ac¬ 
count  of  the  nonlinear  effects. 

Fnmi  Fig.  9,  we  see  that  the  amplitude  of  the  axial  motion 
increases  with  an  increase  of  the  rotor  runout.  In  all  three 
runout  cases,  the  axial  motion  takes  place  about  an  average 
value  of  Z  which  is  always  positive.  Hence,  the  average  seal¬ 
ing  gap  is  always  larger  than  the  design  clearance  C,. 
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Fig.  7(a)— Minimum  film  ttilcknaas,  (fi/Colmn.  vs  rotor  runout,  y,.  for  var¬ 
ious  typos  of  olastic  support. 


Fig.  7(b)— Rolative  tilt,  y,  vs  rotor  runout,  y„  for  various  typos  of  olsstic 
support. 
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7(c)— Axial  diaplaeomont,  Z,  vs  rotor  runout,  y„  for  various  lypss  of 
slaslle  support. 
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Fig.  8(a)^Minimum  film  thickness,  (h  Co)mM.  vs  rotor  runout,  Yr,  and  cone 
height,  He- 
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till'  \<‘4U«]  Y  i.iti  iu'  (  oiiilMix-d  <){  \('<t(ii  7,/  uliuh 

4)11  7  ,  .ill HU'.  .UK I  .1  i  4»i.iitn^  \<*(  i(ii  7  .  u  hu  li  <)<•- 
priuis  <»M  7,  .ilnuf.  In  this  4  .i\c,  ilu-  maxmiuni  value  ol  7,  is 
ilie  Mini  7  ,  i-  7./  .iiul  ihe  minmumi  value  of  7.  is  7,,  -  7,/ 

1  Ills  \i<*l(ls  an  .iveia^e  value  dial  is  eijual  to  7.,  aiui  an 
anijiliuule  dial  is  eijual  to  7,/.  I  lie  amplitudes  (ti  7.  in  Fig. 
Id  4  01  ifspondmg  to  7,  d  d  7;*)  and  I  aii'i  d.  Idd.  d.  if)d. 
aiul  d  ld‘F  lespei  ti\(‘l\ .  I  he  analvtua!  value  of  7./  found 
hv  using  (/)  is  7,/  '  d.lL'd  whuh  is  about  2d  jieueni  less 
dian  (he  iium<-iKal  amplitude.  1  lie  aver.ige  v.ilues  f)J  7,  m 
Fig.  Id  < on esponding  to  7,  d..S,  d.?.”).  and  1..*)  are:  d.S12, 
(F77  l,aud  I  a  id.  respei  tivelv .  I  he  .inalvti(al  values  (tl  7., 
<al4iiiaU‘<i  hv  die  small  peitiiihalion  anaivsis  loi  the  three 
7,  values  are:  ('  S2i.  d.Sd.P  and  I.fitff).  respetlivelv.  Henee. 
die  alia’  .lit  al  7.,  values  ai  t onlv  about  4  pereeni  higher  than 
die  iiiiuieiu.il  .iverage  values  of  7,. 

We  •  an  tout  hide.  dieielt‘ie,  that  F.tj.  1  1 1  is  v.iliti  in  hodi 
die  linear  and  nonline.ii  an.ilvses  .uul  tluit  theie  is  a  good 
(oiit'i.tlion  belvveen  the  two  cases. 

I  he  phase  angle  between  tiie  vet  (01  s  7,  and  7.  (see  Fig.  S) 
is  sliowii  ill  Fig.  I  I  .IS  .1  hint  iit>n  ol  lime.  For  the  large  runout 
values  (7,  “•  d.7.')  and  7,  “  1.5)  the  phase  angle  is  negative 
througlioui  ilie  period,  iliat  is.  tlie  stator  till  vetior.  7..  al¬ 
ways  lags  tlie  rotor  till  vetior.  7,.  For  the  small  nmtiiii, 
7,  -  d.:b  (he  fihase  angle  is  negative  over  most  of  the  period 
hui  here,  altcrnaielv.  tlie  lilt  of  the  rotor  and  that  of  the 
stator  is  in  the  lead.  I'his  behavior  was  tibserved  experi- 
inentallv  and  reported  in  {S). 
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Fig.  a(e)— Axial  dlaplacamant,  Z,  va  rotor  runout,  y„  and  eons  height,  H. 


hr.iiii  I'i);.  III.  it  .iiii  In’  scrii  lliiil  ihc  iiiiiplilii.lc  iil  7,  is 
.ilinosi  ituU’|K'ii<lt'iil  <>l  tlie  rtiii.nil  7,  .vliilc  ihf  iivfiiigt'  viiliic 
.»(  7,  is  siroiiKly  .iKeitetl  by  7,.  Kcfei  riiiK  to  Viff.  .S,  we  .an 
assiinie  that,  similar  to  the  small  |)ertiirbali(in  analysis  (/), 


CONCLUSION 

A  (onipuier  |)H)({i  iim  was  used  lo  solve  the  complete  set 
o(  the  noiiliiKMi  equations  of  motion  of  the  flexibly  mounted 
stator  in  a  noncontat  ting  (oned-face  mechanical  seal.  Rotor 
runout  as  well  as  assenihly  tolerances  in  the  rorm  of  initial 
stator  misalignment  weie  a. counted  lor.  ('iivitalioti  in  the 
thud  lilin  between  the  in.iting  laces  was  also  taken  into  ac¬ 
count.  The  .oinpiitei  progiain  enables  transient  anaivsis  of 
the  seal  dvnamics  and  provides  simulatioti  of  tlie  seal  t)e- 
havior.  Both  stability  tliresliold  ami  stea.iv-slate  tesponse 
of  the  flexibly  mounted  stator  were  investigated. 

In  geneial.  it  was  loimd  that  the  critical  shaft  s|)eed  cor¬ 
responding  to  stability  threshold  is  quite  high.  Hence,  the 
rivnamic  sialrility  .should  not  he  a  problem  in  the  majority 
ol  nomrintacting  seals.  ,'\  tnore  practical  problem  in  non¬ 
contacting  seals  is  the  steady-slate  dynamic  response  of  the 
stator  resulting  frum  rotor  runout  and  assembly  tolerances, 
i'he  stator  stea.ly-staie  response  is  a  nonsynch ronous  track¬ 
ing  of  the  rotor  runout  with  cyclic  variations  in  the  seal 
clearance,  the  relative  misalignment  between  its  mating  laces, 
and  in  the  pha.se  angle  Irelween  stator  and  rotor  tilts.  A 
parametric  investigation  was  performed  lo  explore  the  ef¬ 
fect  of  various  design  parameters  and  operation  conditions 
on  the  seal  behavior.  A  critical  runout  was  found  that,  if 
exceeded,  results  in  seal  failure  due  to  I.Hal  face  contact.  It 
was  also  found  that  seals  may  fail  Itecause  of  an  excessive 
leakage  resulting  from  an  increase  in  the  seal  dearance 
much  iM-fore  the  i  rili.  al  1  iinoul  is  reached.  This  .  an  hap|>en 
ill  seals  with  low  stiffness  of  the  fluid  him  and  the  flexible 
siipijort.  Small  design  clearance,  high  pressure  differential. 
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fig.  9— Time  variation  of  the  axial  displacement,  Z,  at  an  initial  stator 
misalignment,  y««  =  3f  various  rotor  runouts,  y,. 
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Fig.  10— Tim*  variation  ol  tha  nutation,  y„  at  an  Initial  atator  mlaallgn- 
manl,  yu  =  3,  ond  varioua  rotor  runouta,  y,. 


and  o|]limuin  coning  provide  high  stil  lness  of  the  flni<l  liltn 
.ind  rednie  the  danger  ol  a  tiKi  high  o|>eraling  i  learaiue. 

The  results  ol  the  numerical  analysis  were  totnparerl  with 
live  results  ol'  a  small  perturbation  analysis  that  provides 
mtich  simpler  t  losed  lorm  analytical  .solutioti.  Very  good 
correlation  was  lound  lictween  the  two  analyses  for  tnosi 
cases  ol  practical  applications.  Fair  correlation  was  found 
even  in  cases  which  are  clearly  ottt  of  the  range  of  small 
perturbation.  In  these  cases,  the  small  perturbation  analysis 
vields  resitlts  that  are  on  the  safe  side  with  regard  to  the 
cri'it  al  rotor  rtmoul.  If  a  practical  litnit  is  set  for  the  per¬ 
missible  tnitiitntttn  liliti  thi<  kness  at,  say,  2.')  |K’rcetit  of  the 
design  (le,iran(e.  theti  the  an.il\ii<.il  lesults  of  the  small 
peitinbation  ,iti.ilvsis  can  prohahiv  be  valid  over  the  lull 
range  of  desigti  parameters  and  opeiatioti  (onditions. 
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